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Aims. We determine carbon, oxygen, and their isotopic abundances based on CO and OH spectra in 23 red giant stars, and identify 
possible origin of difficulty in abundance analysis of cool luminous stars. 
>— ^ I Methods. We apply the line-by-line analysis based on the classical micro-turbulent model and 1 D model photospheres. 

^^ , Results. We found purely empirically that there is a critical value of log W/v ~ -4.75 (W is the equivalent width and v the wavenum- 

ber) above which the observed lines do not follow the classical line formation theory based on the micro-turbulent model and that the 
classical abundance analysis can be applied only to the weak lines of log W/v S -4.75. The carbon, oxygen, and their isotopic abun- 
dances in 23 K - M giant stars obtained from such weak lines are roughly consistent with the recent evolutionary models, although the 
'^C/'^C puzzle remains unsolved. On the other hand, it is already known that the very strong lines of log W/v > -4.4 are contaminated 
Q ■ by the contribution from extra molecular layers. The less strong but saturated lines of -4.75 < log W/v S -4.4 (intermediate-strength 

f-^ ' lines) cannot be understood at all with the abundance and turbulent velocity based on the weak lines. From the behavior shown by 

C/3 I these lines and considering other observations such as the detection of H2O lines, not only in the late M giants but also in the early 

C^ , M and K giants, we found that the intermediate-strength lines are essentially the same as the strong lines in that they also include 

contamination from the extra molecular layers. Thus the behavior of the intermediate-strength lines, including those with LEP as high 
as 2 eV, appears to be nothing but a manifestation of the warm molecular envelope or MOLsphere. 

Conclusions. The infrared spectra of K - M giant stars are a hybrid of at least two components originating in the photosphere and 
MOLsphere. Only the weak lines, mostly high excitation, are relatively free from the effect of MOLsphere and can be used to extract 
photospheric abundances. The strong lines and the intermediate-strength lines which actually dominate the observed infrared spectra, 
are badly disturbed by contamination from the MOLsphere. For this reason, they provide little information about the photosphere, 
but instead can be new proves of the warm molecular envelope for which little is known yet. In the interpretation and analysis of the 
infrared spectra of K - M giant stars, it is essential to keep their hybrid nature in mind. 






Key words. Line: formation - stars: abundances - stars: late-type - stars: atmospheres - stars: mass-loss - Infrared: stars 

1 . Introduction considering the hydrodynamics of convection in 3 dimension (3 

D model) (e.g. Asplund 2005). In the regime of the coolest stars 
The basic prmciple of stellar abundance determination was es- ^^^^ ^^ ^^^ gj^^^ ^„j supergiant stars, however, the situation may 
tablished already about half a century ago (e.g. Unsold 1955), ^^ ^ ^-^ ^-^^^^^^^ ^^^ ^^^^^ ^^^ ^^ ^^^^^^^. p^^^j^^^ 
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S but it is only recently that high quality spectroscopic data ai-e ~ ' ''^;,~^^ sTai^s "representing tlie evolmion from the red gi- 

made available for accurate abundance analyses. Actual abun- ^^^ ^^^^^^ ^^^^^ ^^ ^^^ asymptotic giant branch (AGB), CNO 

dance determinations have largely been done using the classical ^^^ ^^^.^ ■ -^ abundances experience drastic changes and ac- 

methods based on grids of model photospheres (e.g^ Gustafsson ^^^.^^^ abundance determination should be vital not only to test- 

et a . 1975; Tsuji 1978; Kurucz 1979; Plez et al. 1992) In many . ^j^^ ^^^jj^^. evolutionary models but also to understanding the 

applications local thermodynainical equilibrium (LTE) was as- chemistry of interstellar matter which is enriched by the mass- 

sumed to hold (e.g. Unsold 1955) and favored for its simplicity j^^^ ^^^^ ^^^ .^^^ ^^^^ ^ .^^ ^^^^ importance, few abun- 

in abundance determinations. The next inajor step was to treat ^^^^^ ^^^^ ^^^^^ ^^ ^.^ resolution spectra have been done 

hne formation without the restriction of the LTE assumption ^^^ ^^^j .^^^ ^^^^.^^ compared to other type of stars. It should be 

(e.g. Mihalas 1978) More recently, the classical model photo- rei^embered that high quality infrared spectra needed for this 

spheres describing the structure of the photosphere only in the j^^^^ ^^^^ available, thanks to the Fourier transform 

radial direction (1 D) are being challenged by the new models spectroscopy (FTS) pioneered by P. & J. Connes (e.g. Connes 

1970) and developed further in the 1970's (e.g. Maillard 1974; 

Send ojfprmt requests to: T Tsuji j^^jj ^^ ^j ^^^^. rjj & g^^^j^ ^984)^ ^^^^ before high ef- 

Tables 2 and 3 are only available m electronic ^ • ^ 1 ^ ^ u r^r^r^ ■ » j »• 1 
„ , „_ „ . ■' ^ , ncient detectors such as CCD innovated optical spectroscopy, 
form at the CDS via anonymous ftp to cdsarc.u- r^, , , , , • , ,■,-,, , 
strasbg.fr (130.79.128.5) pub/A+A/489/1271 or via Thus the problem here was not m the quality of the observed 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+ A/489/127 1 The FTS data but should be due to some inherent difficulties in the spec- 
spectra used in this smdy are also available from the same location tral analysis of very cool stars. So far, some efforts to determine 
(subdirectory: ftsdata/), by courtesy of K. H. Hinkle and S. T. Ridgway the CNO and their isotopic abundances were made for oxygen- 
(Kitt Peak National Observatory). rich case (e.g. Harris & Lambert 1984; Smith & Lambert 1985, 
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1986, 1990; Harris et al. 1985; Tsuji 1986, 1991; Aoki & Tsuji 
1997) as well as for carbon-rich case (e.g. Lambert et al. 1984; 
Harris et al.l987; Ohnaka & Tsuji 1996, 1999; Abia & Isern 
1997; Ohnaka et al. 2000; WaMin et al. 2006; Abia et al. 2008). 
Although many interesting results were shown by these works, 
we were not convinced entirely that the stellar abundance analy- 
sis could be done consistently for cool giant stars. 

For example, in our previous analysis of the CO first over- 
tone bands in M giant stars, we found that most lines could 
be interpreted consistently with the classical model in terms of 
the micro-turbulent velocity and carbon abundance (Tsuji 1986, 
hereafter Paper 1). The very strong lines of low excitation could 
not be fitted into this scheme but it was shown that these strong 
lines were contaminated by the contributions of outer molecular 
layers (Tsuji 1988, hereafter Paper II). We then analyzed the CO 
second overtone bands of the same sample of M giant stars and 
found that the resulting carbon abundances appear to be larger 
by about a factor of 2 compared with those based on the CO first 
overtone bands (Tsuji 1991, hereafter Paper III). The origin of 
these discrepant results remains unsolved yet. The difficulty in 
abundance analysis of very cool stars was also noticed by Smith 
& Lambert (1990), who found that the strong OH 1-0 and 2- 
bands gave unreasonably large oxygen abundances compared 
with those based on the weaker bands and that strong CO lines 
cannot be used for abundance analysis as well. 

The carbon and oxygen abundances in other stars are by no 
means well established yet. Even for the sun, the classical val- 
ues widely adopted (e.g. Anders & Grevesse 1989) had to be 
challenged by a new approach based on the 3 D models, and a 
reduction of carbon and oxygen abundances by more than 50% 
from the classical values was suggested (Allende Pietro et al. 
2002). The classical values, however, were defended by a de- 
tailed analysis of CO lines with semi-empirical model solar pho- 
tosphere by Ayres et al. (2006). These results remind us that the 
accurate modeling of the photospheric structure is a critical is- 
sue in abundance determinations. It is unknown if the difficulties 
in abundance analysis in M giant stars outlined above are due 
to a limitation of the classical 1 D model photosphere for cool 
giant stars or if they are related to other problems. For example, 
spectra of red giant stars revealed such an unexpected feature as 
H2O lines, not only in the early M giants and K5 giant a Tau 
(Tsuji 2001) but also in the earlier K giant star a Boo (Ryde et 
al. 2002). These features cannot be fitted into our present un- 
derstanding of the photospheres of red giant stars, and suggest 
a new problem in the photosphere and/or in the atmosphere ex- 
tending beyond the photosphere. For this reason, we analyze the 
nature of the difficulty in interpreting the spectra of cool giant 
stars in some detail in this paper. 



Table 1. Program stars and their basic stellar parameters 



Obj(BS/HD) 


Sp. type 


^ 

-'ctf 


Mj 


Mass '■ 






K 


mag. 


Me 


a Tau ( 1457 ) 


K5+III 


3874 


-1.7 ±0.2 


1.5 ±0.3 


S Oph ( 6056 ) 


M0.5III 


3790 


-2.2 ± 0.3 


1.6 ±0.3 


vVir(4517) 


MIIII 


3812 


-2.2 ± 0.4 


1.7 ±0.4 


Q'Cet(911 ) 


M1.5IIIa 


3909 


-3.2 + 0.3 


3.6 ± 0.4 


0- Lib ( 5603 ) 


M2.5III 


3596 


-3.4 ± 0.5 


2.2 ± 0.5 


A Aqr ( 8698 ) 


M2.5III 


3852 


-3.4 ± 0.7 


3.7 ± 1.2 


/3 Peg (8775) 


M2.5II-III 


3603 


-3.3 ± 0.2 


2.2 ± 0.3 


T* Eri ( 1003 ) 


M3+IIIa 


3712 


-2.9 ± 0.4 


2.0 ± 0.4 


H Gem ( 2286 ) 


M3III 


3643 


-3.3 ±0.3 


2.3 ± 0.5 


5Vir(4910) 


M3III 


3643 


-2.4 ± 0.3 


1.4 ±0.3 


10 Dra ( 5226 ) 


M3.5III 


3730 


-2.9 ± 0.3 


2.1 ±0.4 


p Per (921 ) 


M4II 


3523 


-4.1 ±0.4 


3.2 ±0.5 


BS6861 (6861 ) 


M4 


3600 


-5.2 ± 2.0 


6.3 ± 4.0 


(52Lyr(7139) 


M4II 


3420 


-5.5 ± 0.8 


5.5 ± 2.0 


RR UMi ( 5589 ) 


M4.5III 


3397 


-3.4 ± 0.3 


1.6 ±0.3 


a Her ( 6406 ) 


M5Ib-II 


3293 


-5.8 ± 1.6 


5.0 ±2.0 


OP Her ( 6702 ) 


M5II 


3325 


-4.4 ± 0.8 


2.3 ± 1.0 


XY Lyr ( 7009 ) 


M5n 


3300 


-5.1 ± 1.1 


3.7 ± 1.5 


RLyr(7157) 


M51I1 


3313 


-4.3 ± 0.3 


2.1 ±0.5 


RZ Ari ( 867 ) 


M6III 


3341 


-3.5 ± 0.6 


1.5 ± 0.4 


30g Her (6146) 


M6-III 


3298 


-4.2 ± 0.4 


2.0 ± 0.6 


SWVir (114961) 


M7-III: 


2886 


-4.4 ± 0.9 


1.3 ±0.3 


RX Boo (126327) 


M7.5-8 


2850 


-4.6 ± 0.8 


1.3 ±0.3 



'' based on the infrared flux method. The uncertainty of T^ff for indi- 
vidual object is estimated to be ±100 K. 

*" based on the Hipparcos trigonometric parallaxes and measured 
bolometric fluxes. 

'^ based on the evolutionary tracks of Fig. 1. 

Table 4. CNO abundances used in model photospheres 





case a 


case b 


log Ac- 
log A^" 
logV 


8.35 
8.11 
8.69 


8.08 
8.45 
8.69 



" Tsuji (1991). 

>> Aoki & Tsuji (1997). 

'^ Allende Prieto et al. (2002). 



'^C^'^O, and '^c'^O in the H and K band regions, and those of 
'^OH in the H and L band regions. The logarithms of W/v (W 
is the measured equivalent width and v is the wavenumber of a 
line) are given in Tables 3. 



2. Input Data 

2.1. Observed spectra 

Our sample consists of one K giant and 22 M giant stars listed 
in Table 1 and shown in Fig. 1 . We use high resolution infrared 
spectra observed with the FTS of KPNO (Hall et al. 1979). 
Spectra of the H and K band regions are mostly from the KPNO 
archive, and we also observed some spectra of the H and K band 
regions and most spectra of the L band region. Most of these 
spectra were also used in our previous analyses of CO (Papers 
I, II & III), SiO (Tsuji et al. 1994), and NH & CN (Aoki & 
Tsuji 1997). Some details of the spectra including the resolution, 
S/N ratio, and observing date are summarized in Table 2. We 
measured equivalent widths (EWs) of isolated lines of '^C'^O, 



2.2. Model photospheres 

We apply classical 1 D model photospheres that are essentially 
the same as those adopted in our previous works (Tsuji 1978; 
Paper I - III), except that we now assume that the photosphere is 
spherically symmetric (SS) rather than plane-parallel (PP) and 
some molecular opacity data are updated (Tsuji 2002). We in- 
clude 34 elements in the solar composition (Table 1 of Tsuji 
2002) except for C, N, and O for which we assume two cases 
a and b listed in Table 4 Q. Case a represents the mean C and 
N abundances of early M giants and case b of late M giants 
based on our previous analyses (Paper III; Aoki & Tsuji 1997). 



Abundances are expressed throughout on the scale of log Ah 



12.0. 
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CD 
I 



I 



CM 
I 



O 



CSI 

3.65 




Table 5. Model photospheres 



3.6 



3.55 



3.5 



3.45 



\ogT 



eff 



Fig. 1. Program stars shown by filled circles on the HR di- 
agram and recent evolutionary tracks (Claret 2004) shown by 
solid lines. Stellar masses (in unit of solar mass Mq) are indi- 
cated on some evolutionary tracks. 



The masses, radii, effective temperatures, and log g of our model 
photospheres are listed in Table 5 q We assume masses to be 
one and two solar masses to see the effect of mass. However, we 
used the models with M - 2Mq in the following analyses in 
this paper, since this appears to be closer to the masses of most 
red giant stars we are to study (Table 1). Also, the changes in 
the photospheric structure are rather minor if the masses are in- 
creased beyond 2 Mq, since the differences between the SS and 
PP models tend to be reduced for larger masses. We assumed the 
micro-turbulent velocity in modeling to be 3 km s"' throughout. 
Some characteristics of the resulting models are given in 
Table 5: The surface temperature Tq at tq = 10"^ (tq is the op- 
tical depth defined by the continuous opacity at A - 0.81/im) 
and extension parameter defined hy d = r(To = 10^^)/R. Some 
of our models are compared with those by Plez (1992) in Fig. 2. 
These models are based on the similar assumptions except that 
Plez et al.(1992) applied the opacity sampling method while we 
use the band model method in taking the blanketing effect of 
molecular lines into account. The agreements are generally fair 
apart from at the very surface of the coolest model. In evalu- 
ating the flux, we calculate spectra at a sampling interval of 
0.05 cm ' using a detailed line list including CO (Guelachivili et 
al. 1983; Chackerian & Tipping 1983), OH (Jacquinet-Husson 



^ We refer to a model by its case/mass (in M0)/radius (in Rq) IT^s 
(in K). For example, the first model photosphere in Table 5 is referred 
to as a/1.0/50/4000. 



case 


R 

Rq 


Teff 

K 


logg 


K 


d" 


logF^ 


log< 


M = 


Mq 














a 


50 


4000 


1.04 


2196 


1.06 


12.206 


-2.542 


a 


50 


3900 


1.04 


2185 


1.08 


12.193 


-2.572 


a 


75 


3800 


0.69 


2117 


1.13 


12.185 


-2.609 


a 


75 


3700 


0.69 


2094 


1.13 


12.171 


-2.641 


a 


100 


3600 


0.44 


2023 


1.18 


12.160 


-2.679 


b 


100 


3600 


0.44 


2074 


1.18 


12.159 


-2.677 


b 


100 


3500 


0.44 


2037 


1.18 


12.144 


-2.711 


b 


150 


3400 


0.09 


1857 


1.30 


12.137 


-2.755 


b 


150 


3300 


0.09 


1628 


1.29 


12.118 


-2.787 


b 


200 


3200 


-0.16 


1396 


1.40 


12.102 


-2.825 


b 


200 


3100 


-0.16 


1291 


1.38 


12.064 


-2.842 


b 


250 


3000 


-0.36 


1079 


1.51 


12.031 


-2.866 


b 


250 


2900 


-0.36 


1021 


1.47 


11.974 


-2.868 


b 


300 


2800 


-0.52 


895 


1.58 


11.931 


-2.886 


M = 


2 Mo 














a 


50 


4000 


1.34 


2251 


1.04 


12.200 


-2.536 


a 


50 


3900 


1.34 


2234 


1.04 


12.187 


-2.566 


a 


75 


3800 


0.99 


2172 


1.06 


12.178 


-2.602 


a 


75 


3700 


0.99 


2152 


1.06 


12.163 


-2.634 


a 


100 


3600 


0.74 


2089 


1.08 


12.152 


-2.670 


b 


100 


3600 


0.74 


2142 


1.08 


12.151 


-2.668 


b 


100 


3500 


0.74 


2107 


1.08 


12.135 


-2.702 


b 


150 


3400 


0.39 


2003 


1.13 


12.124 


-2.742 


b 


150 


3300 


0.39 


1884 


1.12 


12.104 


-2.773 


b 


200 


3200 


0.14 


1652 


1.17 


12.083 


-2.806 


b 


200 


3100 


0.14 


1525 


1.16 


12.041 


-2.819 


b 


250 


3000 


-0.06 


1342 


1.20 


11.998 


-2.833 


b 


250 


2900 


-0.06 


1234 


1.19 


11.936 


-2.830 


b 


300 


2800 


-0.22 


1102 


1.22 


11.882 


-2.837 



temperature at Tq = 10"* where Tq is the optical depth defined by 
the continuous opacity ^t A = 0.81/jm. 

d = r(TQ = 10^^)/R is a measure of the extension of the photo- 
sphere. 

Fi is the flux in the L band in unit of erg cm"- s"' per AA = 1 cm. 
Rl = Fhoi/Ft where Fboi is the bolometric flux. 



et al.l999), CN (Cerny et al. 1978; Bausclicher et al. 1988), 
H2O (Partridge & Schwenke 1997), and SiO (Lavas et al. 1981; 
Tipping & Chackerian 1981). Based on the resulting spectra, the 
flux in the L band of the wide-band photometry, Fi, and log Ri, 
where Rl = /^boi/^L with the bolometric flux Fboi = o" T^^ln {cr: 
Stefan-Boltzmann constant), are obtained (Table 5). 

The stellar mass has a significant effects on the surface tem- 
perature Tq and the extension parameter d especially at low val- 
ues of Teff, and we will examine its effects on abundance deter- 
mination (Sects. 4.3 & 5.3). The calculated values of log/?/, dif- 
fer only slightly, however, from our previous values (Tsuji 198 1), 
which would be expected because the L band region is relatively 
free from the heavy line-blanketing effect. For Teff = 3600 K, 
we compare the results for cases a and b, but the differences are 
rather minor (Table 5). 



2.3. Bas\c Stellar parameters 

For our sample of 23 red giant stars, the effective temperatures 
are determined by applying the infrared flux method (Blackwell 
et al. 1980) using the log Rl values from Table 5 and the ob- 
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T 



6000 - 



4000 



2000 




4 



logP, 



9 



Fig. 2. Some of our model photospheres (solid lines) are com- 
pared with those by Plez et al.(1992) (dots). MIMq, T^s, 
logg, and RIRq are as follows. A: 1.0/2750/-0.50/295, B: 
1.0/3200/0.00/166, and C: 2.0/3600/0.50/132. 



served values of log Ri (Tsuji 1981). The resulting values of 
Teff agree with our previous results (Tsuji 1981) in general to 
within 50 K. Our Teft scale is also consistent with the empirical 
ones based on the measured angular diameters by Ridgway et al. 
(1980) and more recent results (e.g. Dyck et al. 1998; Perrin et 
al. 1998; van Belle et al.l999; Mondal & Chandrasekhar 2005). 
We also found that our effective temperatures of individual ob- 
jects agree with the empirical values of these authors, in most 
cases, to within about 100 K; we therefore assume an error bar 
of lOOK for our estimations of T^s- With the Hipparcos paral- 
laxes (ESA 1997lf| and the bolometric luminosities obtained by 
the integration of the SEDs (Tsuji 1981), the absolute bolomet- 
ric magnitudes are determined. The uncertainty in Mboi is dom- 
inated by errors in the parallaxes. The errors in the photometric 
data are estimated to be about ±0.1 mag, including absolute cal- 
ibration errors of a few percent (Bersaneli et al. 1991). 

Our sample is shown on the HR diagram (Fig. 1) and an evo- 
lutionary mass is estimated by comparison with the recent evo- 
lutionary tracks by Claret (2004). Unfortunately, nearly half of 
our sample cannot be accounted for by the models of Claret, and 
we had to extrapolate the evolutionary tracks to estimate stellar 
masses given in Table 1 . The error bars of masses are derived 
from those of Mboi and T^-s. For the objects with T^^ lower than 
~ 3000 K, the extrapolations needed are so significant that the 
uncertainties should be much larger, and we hope that the theo- 
retical evolutionary tracks can be extended to account for these 
objects. 



3. Method of Analysis 

Given a model photosphere characterized by a set of mass, ra- 
dius, and luminosity (or reduced to Teff and log g in PP mod- 
els), a given EW can be interpreted in terms of abundance and 
micro-turbulent velocity, within the framework of the classical 
theory of line formation. Actually abundances of 34 elements 
are assumed in our model photosphere, and what we are to de- 
termine is corrections to the abundances assumed. In our anal- 
ysis, we apply the line-by-line analysis based on the classical 
micro-turbulent model. In this method, we plot the abundance 
correction Alog A obtained for an assumed value of ^micro against 
the log W/v of the line. We expect that larger (smaller) abun- 
dance corrections are needed if smaller (larger) values of the 
micro-turbulent velocity are assumed especially for more satu- 
rated lines, while abundance corrections are independent of the 
assumed ^micro value for very weak lines free from saturation 
effect. The micro-turbulent velocity ^micro is determined so that 
abundance corrections for individual lines show the least depen- 
dence on EWs, (as for detail, see Paper I, Sect. 5.1). 

Also, the probable errors (PEs) of the resulting abundance 
(correction) and micro-turbulent velocity are essentially deter- 
mined from the scatter in the data points plotted on the Alog A - 
log Wjv plane, as explained in detail in Paper I (see its Fig. 2 and 
related explanation). Accordingly, the probable errors are robust 
assessments of the internal errors of our analysis and any result 
that has a large PE should be assumed to have a large uncertainty 
(e.g., some ^micro values with large PEs in Table 6). Apart from 
the internal error, we notice that there should be a systematic er- 
ror due to a selection of the lines. The analysis, which is based 
on the classical micro-turbulent model outlined above, should be 
applied only to lines weaker than a critical value, since otherwise 
a large systematic error will be introduced, as will be shown in 
Sects. 4 & 5. 



4. Reanalysis of CO 

We analyzed the lines of the CO first and second overtone bands 
separately in Papers I and III, respectively, and we summarize 
briefly the difficulties in these analyses. We then reanalyze the 
CO lines with the same observed and spectroscopic data as be- 
fore, but treat the first and second overtone bands together. 



4.1. Previous analysis of CO 

The method outlined in Sect. 3 has been applied to the lines of the 
'^C'^O first overtone bands (see Figs, la-g of Paper I). The anal- 
ysis in Paper I was done using PP model photospheres, but we 
confirmed that the results are essentially the same for SS model 
photospheres of the same T^ff. As an example, we repeat the line- 
by-line analysis of a Her in Fig. 3a: For assumed values of ^micro 
= 3.0, 3.5, and 4.0 kms"', the CO lines of log W/v ^ -4.5 
(those shown by filled circles) show behavior expected from 
the classical micro-turbulent model outlined in Sect. 3, and it is 
possible to find a ^micio value for which the abundance correc- 
tions do not depend on log W/v values. The resulting values of 
^,„i„„ = 3.81 +0.43 km s"' andlogAc = 8.08 +0.27 in fact pro- 
vide a null mean logarithmic abundance correction as confirmed 
inFig.3b0 



' This research has made use of the VizieR catalogue access tool, 
CDS, Strasburg, France, via ADC/NAOJ. 



'^ Some differences with the result for a Her in Table 4 of Paper I 
(^micro = 3.30 + 0.38 kms"' and log Ac = 7.93 ± 0.28), may be due 
to the use of SS instead of PP models in the present analysis. Thus the 
effect of sphericity is certainly not negligible but not very important. 
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Although it appeared to be possible that a self-consistent 
abundance analysis could be done for the lines of log W/v ^ 
-4.5, it was found that the lines of log W/v > -4.5 could not be 
fitted into the above scheme, but behave quite differently (also 
see Figs. 1 and 2 of Paper II). From a detailed analysis of the 
line profiles, strengths, and velocity shifts of these strong lines, 
we suggested that these lines should be contaminated by the con- 
tributions of extra molecular layers above the photosphere as de- 
tailed in Paper II and hence these lines cannot be used for abun- 
dance analysis. 

We then analyzed the lines of the CO second overtone bands 
by the same method as for the CO first overtone bands (Paper 
III). In this case, many weak lines were available and the abun- 
dances determined were relatively insensitive to the value of 
^micro- At the Same time, the value of ^mkro could not be specified 
well for some cases (Fig. 2 of Paper III). The resulting carbon 
abundances based on the CO second overtone bands, however, 
appeared to be about twice as high of those based on the CO 
first overtone bands. The origin of such a systematic difference 
between the carbon abundances based on the CO first and second 
overtone bands is difficult to understand, since both the analyses 
of the first and second overtone bands have been done internally 
self-consistent. 

In view of a general preference of weak lines in abundance 
analysis, it may be possible to adopt the carbon abundances 
based on the CO second overtone bands and abandon the results 
based on the first overtone bands. However, the CO first overtone 
bands are useful for investigating isotopic abundances, since 
'^C'^O and '^C'^O lines of the first overtone bands are good 
indicators of '^C and '^O abundances, respectively (Sect. 6). We 
should then investigate in some detail the problem with the first 
overtone bands if we are to use them for investigating the iso- 
topic abundances. We assumed that the isotopic ratios can be 
determined even if the elemental abundances themselves cannot 
be determined accurately in our preliminary analysis of the iso- 
topic ratios (Tsuji 2007), but we think it necessary to reconsider 
such a viewpoint. 

4.2. The first and second overtone bands 

Although OH analysis to be detailed in Sect. 5 suggests that the 
lines stronger than a critical value of logW/v ^ -4.75 behave 
differently against the weaker lines, we are unable to see this 
effect in the CO lines of Fig. 3, since only a few lines of the CO 
first overtone bands with log W/v ^ -4.75 could be measured 
especially in the later M giant stars (e.g. only 3 lines could be 
measured in a Her shown in Fig. 3). To include more weak lines, 
we tried to analyze together the first and second overtone bands 
of CO. The results for ^micro = 2.4, 3.0 and 3.6 km/s are shown 
in Fig.4a, where the lines of the first and second overtone bands 
of CO are shown by circles and triangles, respectively. All these 
lines are analyzed together, regardless whether they are of the 
first or of the second overtone bands. Unlike Fig. 3a based on 
the CO lines of the first overtone bands alone, it now appears that 
there is a clear distinction between the lines of log W/v > -4.75 
(shown by the open symbols) and those of log W/v ^ -4.75 
(shown by the filled symbols). Only the lines of log W/v 5 -4.75 
show the behavior expected from the classical micro-turbulent 
model of line formation theory noted in Sect. 3. 

It is evident in Fig. 4a that it is impossible to find a tur- 
bulent velocity that provides a consistent abundance correction 
to all the lines measured covering -5.0 < log W/v < -4.2 . 
Thus, the critical value of logW/v w -4.75 found for the CO 
lines is purely empirical as for OH to be discussed in Sect. 5.1. 
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Fig. 3. a) Logarithmic abundance corrections for the lines of the 
CO first overtone bands of a Her plotted against the observed 
log W/v values for assumed values of ^micro - 3.0, 3.5, and 4.0 
kms '. The analysis outlined in Sect. 3 is applied to the lines of 
log W/v < -4.5 (filled circles) with our model b/2.0/150/3200. 
b) Confirmation of the null logarithmic abundance corrections 
forlogAc - 8.08 and^micro - 3.81 kms"', which are the solution 
of the line-by-line analysis of a Her in a). 



Unfortunately, we thought previously (Paper I) that such a crit- 
ical value is logW/v « -4.5 rather than -4.75. In this paper, we 
confirm that this new critical value of logW/v ~ -4.75 should 
apply to all red giant stars that we are to studyo 

From the weak lines which behave as expected for differ- 
ent values of ^mkro in Fig. 4a, we can determine ^micio value 
using the method described in Sect. 3. The solution is found 
to be ^miei-o = 2.82 + 0.23 km s"' and log Ac = 8.40 + 0.13, 
for which the mean abundance correction is null as confirmed 
in Fig. 4b. On the other hand, the intermediate-strength lines 
remain unexplained. We previously concluded that lines with 
logW/v > -4.5, most of which originate in the low excitation 
levels, should include excess absorption originating in the outer 
molecular layers. It is, however, unclear if this interpretation can 



' With the critical value that discriminates the lines well interpreted 
by the classical line formation theory from those that cannot be, we 
distinguish three groups of lines to be referred to as; 1) the weak lines 
(log W/v S -4.75, although these lines are not necessarily very weak 
lines free from saturation effect). 2) the intermediate-strength lines 
(-4.75 < log W/v < -4.4), and 3) the strong lines (log W/v > -4.4, 
these lines are mostly low excitation lines of CO discussed in Paper II). 
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Fig. 4. a) Logarithmic abundance corrections for the lines of the 
CO first and second overtone bands in a Her plotted against the 
observed logW/v values for assumed values of ^micro = 2.4, 3.0, 
and 3.6 kms '. The CO lines of the first and second overtone 
bands are shown by circles and triangles, respectively (model 
photosphere: b/2.0/150/3300). b) Confirmation of the null log- 
arithmic abundance corrections for log Ac = 8.40 and ^micro - 
2.82 kms"', which are the solution of the line-by-line analysis 
of the weak lines (filled symbols) in a). 



be extended to the intermediate-strength lines and we return to 
this subject in Sect. 7.2 after we examine more cases. 

The results of the same analysis on other late M giants show 
similar patterns as shown for 6^ Lyr and 30g Her in Figs. 5a-b. It 
appears that the abundance corrections for the second overtones 
(filled triangles) are often larger than those for the first over- 
tones (filled circles), but we assume that such a difference may 
be within the margins of error of the measurements, including 
the difference in the assumed continuum levels. The case of RX 
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Fig. 5. a) Confirmation of the null logarithmic abundance cor- 
rections for log Ac = 8.22 and ^micro = 2.86 kms"', resulting 
from the line-by-line analysis of the weak lines of CO in 6^ Lyr 
(model: b/2.0/150/3400). b) The same as a) but for log Ac = 8.35 
and ^n,i„o = 2.61 kms"' in 30g Her (model: b/2.0/150/3300). c) 
The same as a) but for log Ac = 8.32 and ^niicro = 3.41 kms"' in 
RX Boo (model: b/2.0/250/2900). 



Boo shown in Fig. 5c is the most difficult case with only 3 weak 
lines, but nevertheless we could find a solution in this case. But 
the PEs both of ^micro and log Ac, estimated by the method noted 
in Sect. 3 are, however, quite large reflecting the poor quality of 
the data available for this latest M giant in our sample. 

We then examine earlier M giants and apply the same anal- 
ysis to the lines of the CO first and second overtone bands of 
/3 Peg, as shown in Fig. 6a. We applied the same method of the 
line-by-line analysis to the first and second overtone bands sep- 
arately in Fig. Ic of Paper I and Fig. 2 of Paper III, respectively. 
The lines of the second overtone bands were all weak lines and 
the results were ^micro - 4. 1+2.8 kms"' and log Ac - 8.24+0.08 
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Fig. 6. a) Logarithmic abundance corrections for the lines of the 
CO first and second overtone bands in p Peg plotted against the 
observed logW/v values for assumed values of ^micro = 1-6, 2.0, 
and 2.4km S-' (model: a/2.0/100/3600), b) Confirmation of the 
null logarithmic abundance corrections for logAc - 8.27 and 
^micro = 2.04kms"', which are the solution of the line-by-Une 
analysis of the weak lines (filled symbols) in a). 
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Fig. 7. a) Confirmation of the null logarithmic abundance cor- 
rections for logAc = 8.38 and ^micro = 2.20kms"', resulting 
from the line-by-line analysis of the weak lines of CO in a Tau 
(model: a/2.0/50/3900), b) The same as a) but for logAc = 8.13 
and ^niici-o = 2.14kms-' in v Vir (model: a/2.0/75/3800), c) The 



same as a) but for logAc = 8.64 and ^„ 
Cet (model: a/2.0/50/3900). 



1.95 km s in a 



(Paper III). Although the lines of the first overtone bands covered 
the range of -5.1 < log W/v < -4.5, all these lines could also be 
analyzed consistently, and the results were found to be ^micro - 
3.12kms"' and logAc = 7.89 (Paper I). Now, exactly the same 
lines of the first and second overtone bands are analyzed together 



and our results for ^„ 



1.6, 2.0 and 2.4 km are shown in 



Fig. 6a. It appears again that no consistent solution can be found 
for all these lines covering -5.3 < log W/v < -4.40. Instead the 



* Note that our previous analysis of the CO first overtone lines with 
-5.1 < logVV/y < -4.5 appeared to provide a solution (Fig. Ic of Paper 
I), but it now appears to be a wrong solution because of the two reasons. 
First, as we have just found, it included the intermediate-strength lines 
to which the line-by-line analysis based on the classical micro-turbulent 
model cannot be applied. Second, it included only few weak lines. With 
such lines of the CO first overtone bands, the line-by-line analysis con- 



lines are divided into two groups separated at logW/v ^ -4.75, 
and a consistent abundance correction can be obtained from only 
the weak lines. The resulting ^nikro = 2.04 + 0.07km s"' and 
logAc - 8.27 + 0.05 are confirmed to provide null abundance 
coiTections as shown in Fig. 6b. The resulting carbon abundance 
agrees well with that of Paper III but differs substantially from 
that of Paper I. It is now clear that the analysis of Paper I in- 
cluded the lines of logW/v > -4.75, which were inappropriate 
for abundance analyses. 



verged to an incorrect solution, and this is an example of a systematic 
error due to an improper selection of lines, as noted in Sect. 3. Such a 
convergence to a wrong solution could be prevented and a correct crit- 
ical value could be found, if a sufficient number of weak lines could 
have been included in our analysis of Paper I. 
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Fig. 9. a) Confirmation of the null logarithmic abundance cor- 
rections for log Ac = 8.41 and ^micro = 3.57 km s"', resulting 
from the line-by-line analysis of the weak lines of CO in r"* Eri 
(model: a/2.0/75/3700), b) The same as a) but for log Ac = 8.26 
and ^n,ic,o = 2.95kms-i in R Lyr (model: b/2.0/150/3300). c) 



The same as a) but for log Ac = 8.26 and ^„ 
SW Vir (model: b/2.0/250/2900). 



: 4.21 km s in 



The results of the same analysis on other early M giants (in- 
cluding a K giant) show similar patterns. Some examples are 
shown for a Tau, v Vir, and a Cet in Figs. 7a-c. It is to be noted 
that the anomalous behavior of the lines stronger than the critical 
value already appears in the K5 giant a Tau, and in all the early 
M giants that we have examined. The nature of the intermediate- 
strength lines, however, is quite difficult to understand. We have 
previously suggested that the lines with log W/v > -4.5 may be 
contaminated by the contributions of the warm molecular lay- 
ers above the photosphere, but this interpretation has so far been 
considered only for low excitation lines observed in late M gi- 
ants (Paper II). It is to be noted, however, that the anomalous be- 
haviors of the intermediate-strength lines shown on the A log A 
- log W/v diagrams are quite similar between the late (Figs. 4 & 
5) and early (Figs. 6 & 7) M giants. 



micro — 

0.07 based on the weak 



Finally, we noticed that the case of p Per is somewhat dif- 
ferent. We again apply the line-by-line analysis for ^micro = 3.0, 
3.5 and 4.0kms"' as shown in Fig. 8a. The resulting ^, 
3.59 + 0.64km S-' and log Ac = 8.27 
lines are confirmed in Fig. 8b. It is to be remembered that exactly 
the same lines of the first overtone bands alone (but omitting 
Unes of logW/v > -4.5) resulted ^micro = 2.68 + 0.30km s"' 
and logAc = 8.00 + 0.22 (see Fig. le of Paper I), although 
it now appears that almost all the lines used in Paper I were 
the intermediate-strength lines, which may be inappropriate for 
abundance analysis. Although the weak lines behave as ex- 
pected, the intermediate-strength lines do not at all in that the 
abundance corrections are now mostly negative, while they were 
mostly positive in the cases examined so far (see Figs. 4 - 7). 
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We found more or less similar patterns in other objects such 
as t'^ Eri, R Lyr, and SW Vir shown in Figs. 9a-c. We note that 
the presence of the warm molecular layers is already known in 
late M giants including p Per, R Lyr, and SW Vir (Paper II) 
and the strong low excitation lines in Figs. 8 & 9 show posi- 
tive abundance corrections consistent with the excess absorption 
in the outer molecular layers, even though weaker intermediate- 
strength lines show negative abundance corrections. Thus, the 
presence of the extra molecular layers in these objects can still 
be considered, but we will return to this subject after we also 
analyze OH lines (Sect. 7.2). 

We conclude that the weak lines of CO behave quite well as 
expected from the classical micro-turbulent model described in 
Sect. 3 and hence photospheric abundances can be determined 
with reasonable accuracy from these lines, although the behav- 
iors of the intermediate-strength lines cannot yet be understood. 
Thus, our initial purpose to determine carbon abundance could 
be achieved, and we summarize the resulting carbon abundances 
based on the lines of the CO first and second overtone bands for 
23 objects in the 6-th column of Table 6. We found that the re- 
sults do not agree at all with those of Paper I but agree rather 
well with those of Paper III except for a few cases. We conclude 
that the carbon abundances of Paper I cannot be correct because 
they were based on the lines that we found to be inappropriate 
for abundance analysis on a purely empirical basis. 

Although our new log Ac values agree rather well with those 
of paper III, we notice for a few objects that there are differences 
as large as 0.3 dex, which have a direct effect on the determi- 
nation of the nitrogen abundances based on CN lines (Aoki & 
Tsuji 1997). Roughly speaking, if Pc and Pn are determined by 
CO and N2 formations, respectively, 

PcN ^ PcPn ^ AcAll\ 

where Pc, Pn, and Pcn are the partial pressures of C, N, and 
CN, respectively. Then, Pcn remains unchanged if a change of 
Ac by AlogAc is compensated for by that of An by Alog An ~ 
-2 AlogAc ■ For a more accurate assessment of the effect of the 
revised Ac on An, we use the line intensity integral fyix) (see 
the Appendix of Paper III) defined so that 

log(W/y)wk = log^/ + r,Oi^). 

We determine a new An, with the revised Ac, so that it produces 
the same Fyix) based on the previous values of Ac (Paper III) 
and An (Aoki & Tsuji 1997). It is impossible to do this analyti- 
cally, but a few trials and errors are sufficient to find a solution. 
We apply this correction in a typical line with the PP models 
used previously for the determination of the N abundances. The 
resulting change in log An is found to be between -Alog Ac and 
-2 AlogAc, and the revised values of log An based on the lines 
of CN Red System are given in the 14-th column of Table 6. 



4.3. Uncertainties 

So far, we have only considered internal errors (PEs) in our re- 
sults shown in Table 6, but the PEs are useful indicators of the 
reliability of the results, as noted in Sect. 3. For example, in- 
spection of Table 6 reveals that the PEs for f micro are larger than 
1 km s ' for BS6861, SW Vir, and RX Boo, and this may be be- 
cause the numbers of weak lines available in these cases are only 
a few. A large systematic error due to an improper selection of 
lines in Paper I is corrected with the recognition of the possi- 
bly true critical value of log W/v = -4.75 in this paper. Other 
possible uncertainties are due to errors in the input parameters 



such as Teff, M, R, and to our modeling of the photospheres. 
We previously examined the effect of uncertainties in T^ff and 
log g within the framework of PP models. The results (Fig. 14 of 
Paper III) indicated that the uncertainties are AlogAc ^0.1 for 
AFeff = 100 K and Alog Ac ~ 0.2 for Alog^ = 0.5. 

The use of SS models, however, introduces another problem: 
Inspection of Table 5 reveals that the photospheric extensions 
and hence surface temperatures show considerable differences 
between models of one and two solar masses especially at low 
Teff. In the model photospheres of Teff - 3000 K, for example, 
the photospheric extension increases from 0.20 R to 0.51 P and 
the surface temperature decreases from 1 342 K to 1 079 K if mass 
is reduced from 2.0 Mq to 1 .0 Mq (Table 5). We assume 2.0 Mq 
in our analysis in this paper throughout, but mass is the most 
uncertain parameter For this reason, we examined the effect of 
mass for a typical case of a Her as an example. We obtained the 
results shown in Table 6 (log Ac = 8.40 and ^micjo = 2.82kms"') 
with a model of 2.0 Mq (b/2. 0/1 5 0/3 3 00). We repeat the same 
analysis with a model of M = 1.0 M© ( b/1. 0/150/3300). The 
resultis log Ac = 8.29+0.13 and^^jcro = 3.05 ±0.24 km s^'. The 
difference in derived carbon abundance is 0.1 1 dex and turbulent 
velocities differ by less than 10 %. Thus, the conclusion obtained 
with PP models remains almost unchanged if we use SS models. 
The reason for this may be due to the use of only the weak lines, 
which are formed rather deep in the photosphere. Also CO is 
quite stable at low temperatures because of its high dissociation 
energy; its abundance may be hardly affected by changes of the 
physical conditions. This is in marked contrast to the case of the 
Sun in which CO is just forming at temperatures considerably 
higher than in M giants and its abundance is highly sensitive to 
temperature. 

4.4. Comparisons witli otiier autlior's results 

We compare our carbon abundances with those by Smith & 
Lambert (1990) for objects analyzed in both analyses presented 
in Table 8. Some significant differences are found for v Vir, 
/3 Peg, and p Per, which show differences of 0.38, 0.15, and 
0.19 dex, respectively. The agreements for other objects, how- 
ever, are generally within the internal errors of both analyses. 
Unfortunately, we are unaware of other high resolution studies 
of photospheric abundances of cool giant stars. 

We note, however, an interesting intermediate-resolution 
analysis for a large sample of 70 M giant stars (Lazaro et al. 
1991). Their result implied rather low carbon abundances consis- 
tent with our previous result of Paper I, which we just discarded. 
Since their intermediate-resolution analysis necessarily included 
the intermediate-strength lines as we also did in our previous 
analysis in Paper I, their result demonstrates that the low and 
high resolution analyses provide similar results if applied to the 
same types of lines. Certainly, intermediate-resolution analysis 
should be useful for applications to faint objects in remote stel- 
lar systems. However, the problem of how to treat the stronger 
lines that have not been well modeled by the classical line for- 
mation theory must be settled before lower resolution spectra 
can be used successfully in abundance analyses. 

Another interesting attempt of medium resolution spec- 
troscopy was reported by Decin et al.(2003), who analyzed in- 
frared spectra observed with ISO SWS. They carefully compared 
the observed and synthetic spectra and estimated basic stellar pa- 
rameters, including the CNO and their isotopic abundances. For 
synthetic spectra of medium resolution, it may again be difficult 
to avoid the use the intermediate-strength and strong lines inap- 
propriate for abundance analysis. Nevertheless the resulting car- 
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Table 6. Carbon and oxygen abundances and micro-turbulent velocities (with probable errors(PEs) in 23 red giant stars 



star 


N " 






IWH* 


log Ac 


Smicro 

kms-' 




IWH* 


N,^ 


J^M,k 


logAo 


Smicro 

kms"' 


log A/ 


a Tau 


13 


13 


25 


9 


8.38 ± 0.04 


2.20± 0.09 


23 


21 


28 


25 


8.79 ± 0.04 
8.76 ± 0.09 


3.85 ±0.13 
3.81 ±0.94 


8.05 


SOph 


5 


5 


21 


6 


8.27 ± 0.08 


2.46±0.16 


42 


27 






8.83 ±0.14 


2.84 ± 0.54 


8.58 


vVir 


6 


6 


22 


10 


8.13 ±0.08 


2.14±0.15 


55 


20 






8.84 ± 0.05 


2.67 ± 0.24 


8.11 


Q-Cet 


6 


6 


11 


3 


8.64 ± 0.05 


1.95±0.08 


38 


18 






8.98 ±0.14 


4.25 ± 1.03 


7.73 


crLib 


12 


12 


14 


2 


8.23 ± 0.03 


2.09±0.13 


51 


19 






8.82 ±0.10 


2.03 ± 0.20 


8.15 


^Aqr 


5 


5 


13 


4 


8.56 ±0.11 


1.97 ±0.15 


30 


7 






8.92 ±0.15 


6.24 ± 6.01 


7.88 


jSPeg 


12 


12 


29 


7 


8.27 ± 0.05 


2.04 ± 0.07 


53 


16 






8.77 ±0.12 


3.12 ±0.72 


8.11 


T^Eri 


6 


5 


20 


1 


8.41 ± 0.04 


3.57 ± 0.46 


25 


9 






8.87 ±0.19 


3.38 ± 1.32 


7.94 


fi Gem 


19 


19 


34 


9 


8.32 ± 0.06 


2.00 ±0.11 


36 


16 






8.81 ±0.13 


2.58 ±0.33 




SYir 


16 


16 


19 


4 


8.50 ± 0.06 


2.20 ±0.11 


36 


11 






8.84 ±0.15 


2.47 ± 0.45 


7.90 


lODra 


5 


5 


15 


3 


8.43 ±0.14 


2.69 ± 0.26 


30 


15 






8.92 ±0.18 


2.44 ± 0.41 


8.11 


pPer 


13 


13 


22 


2 


8.27 ± 0.07 


3.59 ± 0.64 


43 


7 






8.87 ± 0.07 


3.11 ±0.98 


8.15 


BS6861 


5 


4 


12 


2 


8.38 ±0.13 


4.91 ± 1.29 
















S^hyi 


5 


5 


13 


1 


8.22 ±0.18 


2.86 ± 0.40 


19 


3" 










8.41 


RRUMi 


15 


15 


15 


4 


8.26 ±0.12 


2.48 ± 0.23 


34 


5 






8.76 ± 1.05 


2.66 ± 2.60 


8.07 


a Her 


11 


10 


17 


3 


8.40 ±0.13 


2.82 ± 0.23 


22 





42 


16 


8.85 ± 0.05 


2.92 ± 0.09 


8.14 


OP Her 






21 


4 


8.27 ± 0.60 


2.53 ± 0.72 
















XYLyr 






14 


4 


8.15 ±0.15 


2.96 ± 0.24 
















RLyr 


14 


14 


17 


2 


8.26 ± 0.08 


2.95 ± 0.24 


25 


4 


52 


31 


8.61 ± 0.04 

8.62 ± 0.05 


2.99 ±0.10 
6.53 ± 3.20 


8.06 


RZAri 


12 


12 


26 


7 


8.10 ±0.09 


2.54 ±0.15 


52 


21 


28 


15 


8.57 ±0.14 
8.67 ± 0.07 


3.35 ±0.28 
3.49 ± 0.50 


8.44 


30g Her 


6 


5 


24 


3 


8.35 ±0.11 


2.61 ± 0.26 


23 





38 


26 


8.75 ± 0.09 


2.49 ±0.16 


7.92 


SWVir 


11 


6 


15 


1 


8.26 ± 0.20 


4.21 ± 2.38 






29 


16 


8.54 ±0.11 


1.83 ±0.71 


8.46 


RXBoo 


6 


2 


13 


1 


8.32 ± 1.62 


3.41 ± 2.38 






13 


12 


8.55 ±0.10 


2.04 ±0.18 


8.26 



The number of CO lines measured in the H band region, of which the number of weak lines (log W/y S 

The same as " but in the K band region. 

The same as " but for OH lines in the H band region. 

The same as " but for OH lines in the L band region. 

with only 3 weak lines of similar intensities, the micro-turbulent velocity cannot be fixed by our line-by-line analysis 

N abundances based on CN (Aoki & Tsuji 1997) corrected for the effect of the revised C abundances. 



4.75) is given in the next column. 



bon abundances agree rather well with the other results except 
for a Cet, for which our result shows unusually high C abun- 
dance. 



5. The case of OH 

Hydroxyl radical OH shows rich spectra due to ro-vibrational 
transitions in the infrared region. The fundamental and first over- 
tone bands are observed well in the L and H band regions, re- 
spectively, in cool oxygen-rich stars. Spectroscopic data such 
as the line positions and intensities of OH are taken from the 
GEISA databank (Jacquinet-Husson et al. 1999), which covers 
well the high excitation transitions of OH observed in stellar 
spectra. 

5.1. The OH fundamental bands 

First, we applied the same method as used for CO to the lines of 
the OH fundamentals. We include the carbon abundance deter- 
mined in the previous section as an input abundance for each star. 
Since OH formation depends on Aq-Aq rather than on Aq, deter- 
mination of the carbon abundance is a prerequisite for OH anal- 
ysis to determine the oxygen abundance. An example is shown 
for a Her in Fig. 10a, in which the results of the analysis with 
Smicro = 2.0, 3.0, and 4.0 kms"' are shown. We again encoun- 
tered the same difficulty as for the case of the CO first overtone 
bands in that the stronger lines do not follow the pattern expected 



from the weaker lines. Furthermore, it appears very clearly that 
these weaker and stronger lines are divided at log Wjv « -4.75 
and not at log Wjv ~ -4.5 suggested from our previous analysis 
of the CO first overtone bands alone (Paper I). In fact, we first 
noticed that log W/v x -4.75 should be an important value for 
dividing the lines into two groups by this analysis of the OH fun- 
damental bands. With this result in mind, we also confirmed that 
this value is similarly important for the CO lines, but clearly seen 
only if the first and second overtone bands are analyzed together 
(Sect. 4.2). 



As for CO, the weaker lines show the behavior expected from 
the classical micro-turbulent model noted in Sect. 3. However, 
the stronger lines, which should suffer heavier saturation effect, 
do not show the behavior expected for the saturated lines. Thus, 
we apply the same grouping of lines as defined for CO (see foot- 
note 5) to OH. From the weak lines, we found that the abundance 
corrections from different lines show a consistent value for ^mkro 
- 2.92 kms"' as shown in Fig. 10b. We found again that it is 
impossible to smooth out the reflection at log W/v x -4.75 by 
any choice of the micro-turbulent velocity, and thus the critical 
value of log W/v ~ -4.75 can be deemed as well defined em- 
pirically. This result is consistent with Smith & Lambert (1990), 
who found that OH lines of 1-0 and 2-1 bands (mostly stronger 
than the lines of 3-2, 4-3, 5-4 bands) provide unreasonably large 
values of the oxygen abundance and cannot be used for abun- 
dance determinations. 
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0.5 dex. 
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mic ' 
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H 1 1 1 \ 1 \ h 



b) a Her 



2.92 km/s 



S,*o 



_l I L_ 



_L 



-5.5 



-4.5 



log W/u 



Fig. 10. a) Logarithmic abundance corrections for the lines of 
OH fundamental bands in a Her plotted against the observed 
log W/v values for assumed values of ^micro = 2.0, 3.0, and 4.0 
km s-i (model: b/2.0/150/3300). b) Confirmation of the null log- 
arithmic abundance corrections for logAo = 8.85 and ^micro - 
2.92 kms ', which are the solution of the line-by-line analysis 
of the weak lines (filled circles) in a). 
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Fig. 11. a) Confirmation of the null logarithmic abundance cor- 
rections for logAo - 8.79 and ^micro - 3.85 kms"', resulting 
from the line-by-line analysis of the weak lines of OH in a Tau 
(model: a/2.0/50/3900), b) The same as a) but for logAo = 8.61 
and ^micro = 2.99 kms"' in R Lyr (model: b/2.0/150/3300). c) 



The same as a) but for logAo - 8.54 and ^„ 
SW Vir (model: b/2.0/250/2900). 



1.83 kms in 



It is unclear if the same problem exists for the OH lines of 
K giant a Tau, since almost all lines measured are the weak 
lines. We recall that the anomalous behavior of the intermediate- 
strength lines of CO was already noticed in the K giant a Tau 
(Sect. 4.2). We are able to use almost all the OH lines measured 
in a Tau in the abundance analysis simply because they are the 
weak lines and the result is shown in Fig. 11a. However, all the 
other late M giants show similar pattern as for a Her As exam- 
ples, the results for R Lyr and SW Vir are shown in Figs. 1 lb 
and lie, respectively (30g Her, RZ Ari, and RX Boo also show 
similar patterns). 

We assume that the photospheric abundances of oxygen 
could be obtained exclusively from the weak lines, and sum- 
marize the resulting oxygen abundances and turbulent velocities 
obtained from the OH fundamentals in Table 6 for 7 objects for 



which the L band spectra are available (those objects with non- 
zero N'l'' in the 1 1-th column). 

Again, a problem is why the intermediate-strength lines can- 
not be used for abundance analysis. We suggested that some con- 
tributions from the outer molecular layers will be disturbing the 
most strong lines in the case of CO. To examine if this may also 
be the case for OH, we plotted the lower excitation potentials 
(LEPs)of all the lines in Fig. 10a against log W/v in Fig. 12b. For 
comparison, we made a similar plot in Fig. 12a for CO lines of 
a Her analyzed in Fig. 3. Unlike the case of the strong CO lines 
originating from the levels of very low LEP less than 0.5 e V, the 
intermediate-strength lines of OH are not necessarily originat- 
ing from the very low LEP, but rather from the levels with LEP 
~ 0.5-1.5 eV However, the intermediate-strength lines of CO are 
also not originating from the levels with very low excitation lev- 
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log W/v 

Fig. 12. Lower excitation potentials (LEPs) of the lines plotted 
against log W/v values, a) CO first overtone bands observed in a 
Her b) OH fundamental bands observed in a Her. 



els, and the intermediate-strength lines of CO and OH may have 
the same origin as will be discussed in Sect. 7.2. 

5.2. The OH first overtone bands 

Unfortunately, we have only 7 spectra of the L band region, and 
we analyze the OH first overtone bands observed in the H band 
region to extend our sample. Unlike the most molecules that the 
overtone bands are weaker than the fundamental bands, the case 
of OH is rather special in that the first overtone bands are as 
strong as the fundamental bands or even stronger. As a result, it 
is more difficult to measure weak lines in the first overtone than 
in the fundamental bands especially in cooler objects, and we 
cannot analyze the first overtones of OH in some cool objects. 
Fortunately, most of these objects are those for which we have 
the L band spectra and oxygen abundances are akeady deter- 
mined from the OH lines of the fundamental bands (Sect. 5.1). 

We apply the same method as used for the fundamental 
bands and an example is shown in Fig. 13a forp Per For this case 
of the overtone bands, compared to the case of the fundamental 
bands discussed in Sect. 5.1, it is shown more clearly that the OH 
lines are divided into two groups that behave quite differently. 
Even though abundances can be determined rather well from 
the weak lines, it is hopeless to determine abundances from the 
intermediate-strength lines. It is also interesting to note that the 
abundance corrections for the intermediate-strength lines show 
a similar behavior to those for the intermediate-strength lines of 
CO in p Per (Figs. 8a-b). 

In the K giant a Tau, most lines measured are weak lines as 
in the fundamental bands, and we use most of these lines for our 
abundance determination (Fig. 14a). In the M3 giant cr Lib, the 



abundance corrections for the intermediate-strength lines show 
large positive deviations compared with those for the weak lines, 
for which the abundance corrections converge to null for log Ac 
= 8.82 and ^n,i„o = 2.03 kms"' (Fig. 14b). The M5 giant R 
Lyr shows similar pattern as for p Per (Fig. 14c). The resulting 
micro-turbulent velocity is unreasonably large, but the large PE 
implies that such a result cannot be realistic. This is because 
small errors in the measured equivalent widths will give seri- 
ous effect on the determination of ^micro value when only a few 
weak lines are available, and this is a difficulty of our line-by-line 
analysis. Nevertheless, the abundance can be determined rather 
well because of the use of the weak lines, which depend little on 
^micro valuc. The difficulty of measuring weak lines in the first 
overtones of OH is more serious for some late M giants includ- 
ing 6^ Lyr, 30g Her, a Her, SW Vir, RX Boo, and we had to give 
up to analyze the H band spectra of these stars. 

We confirm both in the CO and OH lines that the behav- 
iors of the intermediate-strength lines are quite complicated and 
show different patterns for different objects. For example, the 
intermediate-strength lines, both of OH and CO, do not neces- 
sarily show positive abundance corrections compared with those 
for the weak lines in p Per and R Lyr. It is interesting to note that 
CO and OH lines show similar patterns on the A log A - log W/v 
diagrams if these OH and CO lines are measured from the same 
spectrum as are the cases of p Per and R LyrQ 

We summarize the resulting oxygen abundances based on the 
weak lines of the OH first overtone bands for 15 objects in Table 
6 (those with non-zero A^^'' in the 9-th column). The oxygen 
abundances could be determined independently from the fun- 
damental and first overtone bands for three objects, a Tau, R 
Lyr, and RZ Ari. The resulting oxygen abundances, compared 
directly in the 12-th column of Table 6, agree well within the in- 
ternal errors of both the analyses. The micro-turbulent velocities 
also agree well except for R Lyr, for which the reason for the dif- 
ficulty was discussed already (see Fig. 14c). The micro-turbulent 
velocities of a Tau based on the fundamental and overtone bands 
agree quite well (^mkro ~ 3.8 km s"'), but disagree with the re- 
sult obtained from the CO lines (^micro ~ 2.2 km s"'). This result 
reminds us a question what the micro-turbulent velocity means. 



5.3. Uncertainties and comparison witti otiier autliors 

We notice in Table 6 that the values of micro-turbulent velocity 
are sometimes unexpectedly large both for CO and OH, and that 
most cases of ^mid-o ^ 4kms"' are found when the numbers of 
weak lines used are below 10. The ^miao value by our analysis 
based on the method outlined in Sect. 3 is susceptible to small 
errors in the measured EWs especially if only a small number 
of lines can be used. The poor quality of the data is reflected 
in the large PE of the derived ^micio in Table 6, and the resulting 
^micio value cannot be reliable for such a case. For example, ^miao 
values as large as 6 km s ' for A Aqr and R Lyr are almost useless 
because of the very large PEs. 

Unlike CO, OH may be more susceptible to a change in the 
physical condition. To compare with the case of CO (Sect. 4.3), 
we again examined the case of the cool giant a Her. We repeat 
the same analysis of the OH lines of a Her done with a model 



' While the intermediate-strength lines of CO in R Lyr (Fig. 9b) 
and SW Vir (Fig. 9c) show negative abundance corrections, the 
intermediate- strength lines of OH fundamentals in R Lyr (Fig. 1 lb) and 
SW Vir (Fig. 1 Ic) do not. But OH fundamental bands were not observed 
simultaneously with CO bands (see Table 2) and it is possible that line 
formations took place at different conditions in these variable stars. 
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Fig. 13. a) Logarithmic abundance corrections for the lines of 
the OH first overtone bands in p Per plotted against the observed 
log W/v values for assumed values of ^micro = 2.0, 3.0, and 4.0 
km s-i (model: b/2.0/100/3500). b) Confirmation of the null log- 
arithmic abundance corrections for logAo = 8.87 and ^mici-o - 
3.11 kms ', which are the solution of the line-by-line analysis 
of the weak lines (filled circles) in a). 
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Fig. 14. a) Confirmation of the null logarithmic abundance cor- 
rections for logAo - 8.76 and ^micro - 3.81 kms"', resulting 
from the line-by-line analysis of the weak lines of OH in a Tau 
(model: a/2.0/50/3900), b) The same as a) but for logAo = 8.82 
and ^micro = 2.03 kms"' in cr Lib (model: a/2.0/100/3600), c) 
The same as a) but for logAo - 8.62 and ^micro - 6.53 km"' in 
R Lyr (model: b/2.0/150/3300). 



of M = 2.0 M© (b/2.0/150/3300) giving the results shown in 
Table 6 (logAo = 8.85 and ^n,i„o = 2.92kms"'). Now, with a 
model of M = 1 .0 Mq (b/1 .0/150/3300), the results are log Aq = 
8.91 + 0.13 and ^n,i„o = 3.30 + 0.08 kms"'. The differences are 
again unexpectedly small. The reason for this may again be due 
to the use of only the weak lines, which are formed rather deep 
in the photosphere. 

Finally, we again compare our results with those by Smith & 
Lambert (1990) and by Decin et al. (2003) in Table 8. For most 
cases of oxygen abundances, the results agree reasonably well 
within the error bars of these analyses. It is interesting that the 
agreements are reasonably good despite possible differences in 
the model photospheres used. 



6. Carbon and oxygen isotopic ratios 

6. 1 . Line-by-line analysis 

We applied the line-by-line analysis to the isolated lines of 
'2C"'0, '^C'^O, and '^C'^'O. We use only the weak lines but 
include '^C'^O lines of the first and second overtone bands to- 
gether The ^micro valucs determined from the '^C'^O lines of the 
first and second overtone bands in Sect. 4.2 are used to determine 
abundances of '^C'^O and '^C'^O as well. Given a ^mkro value, 
abundance can be determined directly from each line. Some ex- 
amples of the results are shown in Fig. 15a-d for 10 Dra, RZ Ari, 
6 Oph, and jj Gem. Unfortunately, the number of lines of the iso- 
topic species measured is rather small and the weak lines could 
not be measured at all for some cases especially for '^C'^O. 
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Table 7. Isotopic ratios (with PEs) in 23 red giant stars 



Obj. 


WwkdS)" 


A'wkCn)* 


'^C/'^C 


160/170 


a Tau 


2 





10.6 ± 1.0 


>1000 


^Oph 


3 


2 


11.1 ±0.9 


387 ± 68 


yVir 


5 





8.7 ±1.3 


>2000 


aCet 


2 


2 


11.1 ±0.8 


586 ± 47 


cr Lib 


2 





7.5 ±0.3 


>1500 


AAqi 


2 





7.9 ± 1.4 


>1000 


jSPeg 


3 





7.7 ± 0.5 


>2500 


T^Eri 


1 


1 


12.4 ± 0.3 


687 ± 14 


fi Gem 


6 


1 


10.5 ± 1.2 


798 ± 73 


5Vir 


4 





12.3 ± 1.2 


>2500 


lODra 


3 


3 


14.8 ± 1.6 


151 ± 11 


pPer 


2 





9.7 ± 1.0 


>1000 


BS6861 


1 





48.5 ± 2.9 


>1000 


(5-Lyr 


2 


2 


16.2 ± 1.5 


465 ± 41 


RRUMi 


1 





10.0 ± 0.8 


>2000 


a Her 


2 


4 


11.1 ±0.7 


102 ± 8 


OP Her 


2 


1 


11.3 ± 1.2 


329 ±31 


XYLyr 


1 


3 


15.0 ± 0.4 


223 ± 16 


RLyr 


1 


3 


6.4 ± 0.3 


368 ± 44 


RZAri 


3 


2 


7.9 ±0.8 


607 ± 48 


30g Her 


2 


3 


12.5 ± 1.1 


211 ±42 


SWVir 


2 


2 


22.0 ± 4.7 


432 ± 37 


RXBoo 





3 




233 ± 20 



" number of weak lines of '"'C'^O. 
'' number of weak lines of '^C' 'O. 



For example, the lines measured for 10 Dra are '^C'^O (2, 0) 
R25, R27, and R 32 on the spectrum shown in Fig. 16a. This is a 
favorable case in that three '^C'^O lines could be measured, and 
these lines are used for the analysis of Fig. 15a. In the case of RZ 
Ari shown in Fig. 16b, the '^C'^0 lines are weaker and we could 
measure only two lines (R 25 and R 32) for the line-by-line anal- 
ysis of Fig 15b. The '^C'^O lines are still weaker in A Aqr and 
RR UMi shown in Figs .16c and 16d, respectively, and no line 
could be measured in these cases. The difficulty of measuring 
weak lines of '•'C'^O is more or less the same, although at least 
one line could be measured in most objects except for RX Boo. 

The '^C/'^C and '*0/'^0 ratios based on the resulting 
i2q16q^ i3(^16q^ ^jj^ iiQiiQ abundances are given 4- and 5- 
th columns, respectively, of Table 7 with their probable errors 
(PEs). The resulting '^C/'^C ratios are mostly larger compared 
with our previous result (Tsuji 2007), and this is because '^C'^O 
abundances are increased due to the use of the lines of the first 
and second overtones together 



6.2. Synthetic spectra 

We apply the spectral synthesis method as a check of the re- 
sults given in Sect. 6.1. With the carbon abundance and micro- 
turbulent velocity from Table 6 and with the isotopic ratios from 
Table 7, we generate synthetic spectra between 4279 and 4294 
cm"', where '^C'^O (2-0) 7? 25 - 33 lines are found. Before we 
compare the synthetic spectra with observed ones, we corrected 
for the effects of the instrument broadening with the Norton- 
Beer's apodization function 2 (Norton & Beer 1976), which was 
used in apodization of the observed spectra (see Paper I), and 
also of the macro-turbulent broadening based on the measured 
intrinsic line widths (see Table 4 of Paper Illfl As an exam- 



ple, we compare the observed spectrum of 10 Dra (dots) with 
the synthetic one (heavy solid line) in Fig. 16a. In this region, 
the weak Une of '^ci*-© (2-0) /?92 can clearly be seen at 4292 
cm ', which serves as a check of the carbon abundance 0. The 
lines of i^c'^o (2-0) R 25, 27, and 32 used to determine '^c'^o 
abundance in Sect. 6.1 can be reproduced well by the synthetic 
spectrum. Some other lines of '^C'^O such as (2-0) R 26 and 29 
should include extra components, since the synthetic spectrum 
computed without '^C'^O at all (light solid line) shows some 
absorption features. Also some '^C'^O lines are badly blended 
with other lines, but all the lines of '^C'^O are not contradicting 
with the synthetic spectrum computed for '^O/'^O ratio of 151 
(Table 7). As another example, the case of RZ Ari is shown in 
Fig. 16b. Although the '^C'^O lines are weaker than in the case 
of 10 Dra, some unblended features such as (2-0) R25 and 32 
and all the other lines of '^C'^O are not contradicting with the 
synthetic spectrum computed for '^O/'^O ratio of 607 (Table 7). 
The major purpose of computing the synthetic spectra, how- 
ever, is to estimate the '^O/'^O ratio for the cases in which the 
'^C'^O lines are too weak for EWs to be measured accurately. 
For example, all '^C'^O features are too weak for their EWs to 
be measured in the case of A Aqr shown in Fig. 16c and these 
features, including relatively undisturbed ones such as (2-0) R 32 
and 33, do not contradict with the synthetic spectrum computed 
for ""O/'^O ratio of 1000. As another example, the observed 
spectrum of RR UMi shown in Fig. 16d may not contradict with 
the synthetic spectrum computed for '^O/'^O ratio of 2000. This 
conclusion depends largely on subtle depressions in (2-0) R 25 
and 32. On the other hand, definite absorption seen at the posi- 
tion of R 26 and 29 also appear in the synthetic spectrum com- 
puted without '^C'^0 at all (light solid line) and are not due 
to '^C'''0. However, it is possible that the features assumed to 
be due to '^C'^O may still include unknown impurities, and the 
'^O/'^O ratio estimated on such subtle features should at best be 
regarded as a lower estimate to the '^O/'^O ratio. The '^O/'^O 
ratios estimated this way are given in the 5-th column of Table 7 
without error bars. 



6.3. Uncertainties and Comparisons witli other authors 

The serious uncertainties in the analysis of such faint lines as 
'^C'^O in cool stars are the location of the continuum level and 
effect of blending by many weak lines. As a guide to examine 
such problems, we computed synthetic spectra such as Fig. 16 
for the regions with rich lines of '^C'^O (4264 - 4295 cm ') as 
well as of '^C'^O (4235 - 4266cm"') for all the spectra ana- 
lyzed. This was of some help to reject blending lines in our line 
list. Further, we examined the correlations among line intensi- 
ties, line-widths, line-depths, and radial velocities, and rejected 
those lines that show large deviations. The continuum level was 
estimated by connecting the highest peaks by a smooth curve and 
the consistency of the resulting continuum can also be checked 
with the synthetic spectra. However, we must look for the spec- 
tra with still higher resolution to have final check of the adopted 
continuum. Also, a comparison of the different authors' results 
(see below) suggests that the errors in the final isotopic ratios 



** The macro-turbulent velocity is derived by subtracting the effect of 
the micro-turbulent velocity from the line width on the assumption of 



Gaussian micro- and macro-turbulence. The resulting values are; ^,„iicro 
= 2.7, 2.1, 3.9, and 2.8 km s"' for 10 Dra, RZ Ari, A Aqr, and RR UMi, 
respectively. 

' Strong lines noted by the open triangle in Fig. 16d are mostly origi- 
nating from '^C'*0 (3,1) and they are too strong to be used as indicator 
of the carbon abundance. 
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Fig. 15. The ordinate scale shows abundances of '^C'^O ( first and second overtones by filled circles and triangles, respectively), 
'^C'^O (filled circles), and '^C'^O (filled circles) derived from the observed values of log W/v in the abscissa, a) 10 Dra. Note that 
there are three lines of '^C'^O and '^C'^O, of which two are overlapping in the plot, b) RZ Ari. c) S Oph. d) /i Gem. 



due to such uncertainties may be as large as 100% (or factor of 
two). 



The carbon isotopic ratios were studied relatively well and 
we compare our '"C/'-'C ratios with those by Smith & Lambert 
(1990) and by other authors in Table 8. In addition to those cited 
in Table 8, Wallerstein & Morell (1994) reported the '^c/'^c ra- 
tios for 8 early M giants to be 7 - 12, although no object common 
to our sample is included. Our results on '^C/'^C ratios appear 
to agree rather well with the results by other authors in general. 



The determination of the oxygen isotopic ratios is more dif- 
ficult but a detailed identification of '^C'^O features was done 
on the high resolution FTS spectra by Maillard (1974), who es- 
timated '^O/'^O ratio in a Her for the first time. We compare 
our '^O/'^O ratios with those by Harris & Lambert (1984) and 
by other authors in Table 8. Also, Wallerstein & Morell (1994) 
reported lower limits to '^O/'^O ratio in 8 early M giant stars to 
be 85 - 295. Inspection of Table 8 reveals that "'O/'^O ratios by 
different authors differ by factors of 2 - 4. This may largely be 
due to difficulty to locating the continuum levels and to measur- 
ing weak features, and there should be a room to improve the 
results by better spectroscopic data. 



7. Discussion 

7.1. Abundance analysis and photospheric structure 

The abundance analysis of very cool stars is associated with 
many difficulties. In particular, absorption lines tend to be strong 
due to the low continuous opacity at low temperatures of the 
photospheres of cool giant stars. The consequences are severe 
blending and uncertainties in the continuum level, which are 
not only the obvious difficulties by themselves but also make 
it difficult to measure other weak lines which are most impor- 
tant in abundance determination. In fact, if a sufficient number 
of weak lines which are linearly proportional to the effective col- 
umn density of a species under consideration can be measured, 
abundance determination is quite easy. Such an ideal case is sel- 
dom realized especially in cool stars. However, it was thought 
possible to carry out an abundance analysis with the use of the 
micro-turbulent model (or curve-of-growth method) even if sat- 
urated lines had to be used. This model has in fact been applied 
with reasonable success to abundance analysis in general until 
recently. 

The real difficulty in the abundance analysis of cool stars is 
that the classical method based on the micro-turbulent model can 
no longer be applied for lines of equivalent width larger than a 
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Fig. 16. Observed spectra are shown by dots. Some '^C'^O and '^C'^O features are indicated by open and filled triangles, respec- 
tively. Synthetic spectra with assumed ""O/'^O ratios are shown by heavy solid lines while those with null '^O by light solid lines. 
a) Observed spectrum of 10 Dra compared with the synthetic one with '^O/'^O = 151. b) Observed spectrum of RZ Ari compared 
with the synthetic one with "'O/'^O = 607. c) Observed spectrum of A Aqr compared with the synthetic one with '*0/'^0 - 1000. 
d) Observed spectrum of RR UMi compared with the synthetic one with '*'0/'^0 = 2000. 



certain critical value, which we found to be log W/v * -4.75 
for the case of M giant stars (Sects. 4 & 5). Since stronger lines 
are formed in the upper photospheric layers where modeling 
may be increasingly uncertain, this difficulty can be related to 
the problem of the model photospheres. Previously, Smith & 
Lambert (1990) noted that the CO lines stronger than about log 
W/v X -4.8 ~ -4.7 should not be used for abundance analy- 
sis, since these lines are formed in the layers above the contin- 
uous optical depth of t^. ^ 10"^, where the photospheric struc- 
ture cannot be described well by the usual model photospheres. 
We arrive at a similar conclusion that the lines stronger than 
log W/v x -4.75 cannot be used for abundance analysis, but 
by somewhat different argument based on purely empirical evi- 
dence shown most clearly for the OH lines of the fundamental 
bands (Fig. 10). 

Then, a problem is again why the intermediate-strength lines 
show such abundance corrections as large as 1.0 dex or more 
when the weak lines show null corrections? Probably, it is un- 
likely that the uncertainty in the thermal structure of the upper 



layers could produce such a large effect on the intermediate- 
strength lines, since differences in the model photospheres have 
only modest effects on abundances as we have seen before 
(Sects. 4.3 & 5.3). Also, no model of the upper photosphere 
to explain the intermediate-strength lines or the stronger lines 
is known. On the other hand, dynamical structure represented 
by stellar turbulence may have more serious effect on the 
intermediate-strength lines heavily saturated. Taking advantage 
that the line positions and line profiles accurately can be mea- 
sured with FTS, we showed that the line asymmetries and dif- 
ferential line shifts of CO lines are quite large and that a simple 
model of the depth-independent and isotropic Gaussian micro- 
turbulence cannot be justified (Paper 111). 

Even if we know that a simple model of micro-turbulence 
is not sufficiently realistic enough, yet it is unknown how the 
behavior of the intermediate-strength lines (Figs. 4 - 11 & 13 
- 14) can be reproduced with better treatment of stellar turbu- 
lence. In this connection, it is interesting if a more consistent 
treatment of the dynamical structure of stellar photosphere by 
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recent 3 D models can resolve the present difficulty. Recently, a 
detailed numerical simulation of surface convection in red giant 
stars is reported by Collet et al. (2007). Although their simula- 
tions were not extended to the regime of M giants, their result for 
the coolest model of Teft - 4697 K and [Fe/H] - 0.0 indicated 
that a difference between 3D and ID oxygen abundances derived 
from OH lines (not ro-vibrational lines but of A - X system at 
3 150 A) is about 0.1 dex. This seems to be due to the difference 
of the thermal structures between the 3D and ID models, but 
possible effect of the dynamical structure of the 3D models on 
the saturated lines is unknown yet. Their results showed that the 
effect of the 3D structure should be more important in stars of 
the lower metallicities, but it is unclear if such a large difference 
in abundances derived from the weak and intermediate-strength 
lines in our population I sample can be relaxed by considering 
the dynamical effect of convection in the 3 D models. Thus we 
are inclined to conclude that the difficulty we have encountered 
may not be resolved within the framework of the photospheric 
structure. 



7.2. Abundance analysis and extra molecular envelope 

We already know that the strong lines of CO in late M giant 
stars are badly disturbed by the contributions of the extra molec- 
ular layers (Paper II). Now, since the intermediate-strength lines 
and the strong lines behave similarly on the A log A - log W/v 
diagram (e.g. Figs. 4 - 5), an interesting question arises if the 
intermediate-strength lines can also be understood similarly as 
the strong lines. For the late M giants, the intermediate-strength 
lines are not necessarily originating from the levels of low LEP 
(see Fig. 12a), but still increasingly reduced contributions from 
higher excitation lines can be possible for modest kinetic tem- 
peratures of the warm molecular layers, and the trend of the 
intermediate-strength lines together with the strong lines in these 
giants (Figs. 4-5) can be explained consistently. In fact, given 
that the presence of the MOLsphere has been demonstrated in 
the late M giants (Paper II), it should at the same time produce 
excess absorption in the intermediate-strength lines as it does for 



the strong lines. It is then quite natural to assume the same origin 
for the intermediate-strength lines and the strong lines in the late 
M giants. 

Interestingly, the intermediate-strength lines of CO in the 
early M giants (including a K giant) mostly give larger abun- 
dance coiTections compared with those from the weak lines (e.g. 
/? Peg, a Tau, v Vir, a Cet in Figs. 6 - 7). So far, we have not 
considered the effects of outer molecular layers in these early 
M giant staro- But why not? Especially, the behaviors of the 
intermediate-strength lines on the A log A - log W/v diagrams are 
so similar for the late M (Figs. 4-5) and early M (Figs. 6-7) gi- 
ants, it is natural to assume the same origin for the intermediate- 
strength lines in the early and late M giants. Furthermore, we 
identified H2O bands on the ISO spectra of the K - early M gi- 
ant stars included in our present sample (e.g. a Tau, /3 Peg, and 
a Cet) and one possibility is that these H2O bands are originat- 
ing in the outer molecular layers as for the late M giants (Tsuji 
2001). But this is not a unique interpretation: For example, Ryde 
et al. (2002) detected H2O lines on the high resolution spectrum 
of the K giant a Boo and suggested that its origin can be due to 
anomalous structure of the photosphere. We would like to point 
out, however, that the anomalous behavior of the intermediate- 
strength lines and the presence of H2O bands in the K - early 
M giant stars can consistently be understood by assuming the 
presence of the outer molecular layers. 

In the case of CO lines in late M giants such as p Per, R Lyr, 
and SW Vir, and also in the M3 giant t'^ Eri, the intermediate- 
strength lines and the strong lines show quite different behaviors 
as shown in Figs. 8-9: The strong lines show positive abun- 
dance corrections while the intermediate-strength lines show 
mostly negative coiTections for the carbon abundances obtained 
from the weak lines. We first thought that this difference could 
be a reason why we should consider different explanations for 



'" We noticed that the lines of log W/v > -4.5 show excess absorption 
in some early M giants in Paper 11, but we regret that we did not pursue 
this case in more detail. With the critical value of log W/v » -4.5, the 
lines showing excess absorption were limited to only few lines in the 
early M giant stars. 
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the intermediate-strength lines and the strong lines. However, 
we soon arrive at an opposite conclusion that this fact can be 
supporting evidence for the same origin of the intermediate- 
strength lines and the strong lines. In fact, a possible contribu- 
tion from the outer molecular layers can be either absorption or 
emission depending on the extension, temperature, and optical 
thickness of the outer molecular layers for a particular line. We 
can therefore explain the negative abundance corrections for the 
intermediate-strength lines by possible reductions of their EWs 
due to filling in by the emission from the outer molecular lay- 
ers. On the other hand, we have presently no explanation for the 
intermediate-strength lines if we assume photospheric origin for 
them. 

Finally, a problem is whether the intermediate-strength lines 
of OH can be understood in the same way as for the CO lines. 
Unlike CO, OH may not be so abundant in the outer molecular 
layers, since oxygen may be mostly in H2O at relatively cooler 
environment expected for the MOLsphere. For this reason, we 
first thought that the anomalous strengths of the intermediate- 
strength lines of OH may not be due to contamination from the 
outer molecular envelope, but may instead be originating in the 
photosphere itself. However, the infrared OH lines shows be- 
haviors quite similar to the CO lines (e.g. Fig. 4 & 10), and 
it is again most natural to assume the same interpretation for 
the OH and CO lines. Given that the presence of OH is known 
even in the cooler circumstellar envelopes of red giant stars by 
OH masers, it may be possible that OH can be abundant in the 
warmer MOLsphere. Furthermore the assumption of the chemi- 
cal equilibrium should probably not be applied to the outer lay- 
ers and we should rather use the OH lines as probes of non- 
equilibrium processes in the outer layers. Then, we conclude 
that the intermediate-strength lines of OH observed in the late 
M giants such as a Her, R Lyr, and SW vir (Fig. 10 - 11) as well 
as in the earlier M giant such as cr Lib (Fig. 14b) also include 
contributions from the outer molecular layers. 

The possible presence of the outer molecular layers in red 
giant stars is supported by direct measurements of the angu- 
lar extensions. For example, angular diameters of normal K and 
M giant stars, including several objects common to our sample, 
were measured in the strong TiO bands at 712 nm and in the 
nearby pseudo-continuum at 754 nm (Quirrenbach et al. 1993). 
The observed 712nm/754nm diameter ratios appeared to be 
larger than predicted by model photospheres, suggesting that 
the photospheric models cannot describe adequately the outer 
layers of these red giant stars, and the result may lend support 
for the presence of the extra molecular envelope (Quirrenbach 
2001). Also, recent interferometric observations in the infrared 
bands provided direct support for the presence of extended warm 
molecular gas layers not only in Mira variables but also in the 
late M giants (e.g. Mennesson et al. 2002; Perrin et al. 2004). 
On the other hand, the apparent angular diameter of a Her in the 
mid infrared was found to be larger by 30% than its near infrared 
size, but the high resolution 1 1 yum spectrum does not show any 
substantial spectral lines in a Her (Weiner et al. 2003). This re- 
sult was explained by a balance of absorption and emission in a 
spherically extended envelope, and the mid infrared spectra and 
visibilities could be consistent with the warm molecular enve- 
lope model (Ohnaka 2004). 

It is to be remembered that a possible presence of the warm 
molecular envelope in the early M giant stars was previously 
noticed when we first identified H2O bands on the spectrum of 
the M2 giant yS Peg observed with ISO SWS (Tsuji et al. 1997). 
Also, H2O bands were found in several M giants earlier than 
M6 by Matsuura et al.(1999) with IRTS (Infrared Telescope in 



Space) launched in March 1995 by IS AS. These results were 
further extended to several K - M giant stars (Tsuji 2001) and 
also confirmed by other authors (e.g. Decin et al. 2003). The im- 
plication of this finding, however, was not necessarily realized 
very well even by ourselves, and it took some time before we 
recognize that this finding should also be related to the difficulty 
encountered in the abundance analysis of red giant stars. Since 
the temperatures and radial velocities of the MOLsphere differ 
only slightly from those of the upper photosphere, it was diffi- 
cult to recognize the presence of the MOLsphere by the spec- 
troscopic methoco However, we finally recognize the effect of 
the overlaying warm molecular layers on the high resolution in- 
frared spectra of dozens of red giant stars. 

Considering several lines of reasonings outlined above, we 
conclude at last that the presence of a rather warm and quasi 
static molecular envelope above the photosphere is a common 
feature in K and M giant stars. It is quite unexpected and even 
surprising if both the CO and OH lines of intermediate-strength 
with LEP as high as 2.0 eV (Fig. 12) include the absorption 
and/or emission originating in the outer molecular layers. We 
have previously noticed such a possibility mainly for the late M 
giants (Paper 11), but now we find that the presence of the ex- 
tra molecular envelope is a general phenomenon not only in the 
late M giants but also in the early M and K giants as well. Also, 
the fact that the rather high excitation lines are included in these 
intermediate-strength lines implies that the extra molecular en- 
velope is indeed "warm" and should be even warmer than we 
imagined before from the analysis of the strong lines (Paper II). 
Yet we know little about this new feature consisting of warm 
gaseous molecules such as CO, OH, and H2O. We hope it our 
next step to extend our analysis as done for the strong lines 
(Paper II) to the intermediate-strength lines and to clarify the 
nature of the warm molecular layers in more detail. 

We conclude that the observed infrared spectra should actu- 
ally be a hybrid of multiple components originating in the pho- 
tosphere and overlaying warm molecular layers or MOLsphere. 
The hybrid nature of the infrared spectra of red giant stars is now 
confirmed both in the low resolution ISO spectra and high res- 
olution FTS spectra. We must then abandon the traditional idea 
that the stellar spectra are primary defined by the stellar photo- 
spheric structure, at least for the infrared spectra of K - M giant 
stars. 



7.3. CNO abundances 

Even if the photosphere of a cool giant star is veiled by the warm 
molecular envelope, the weak lines of high excitation may re- 
main almost undisturbecQ In fact the weak lines are clearly dis- 
tinguished from the intermediate-strength Unes on the A log A 
- logW/v diagram (Figs. 4 - 11, 13 - 14), and we confirmed 
that the weak lines well behave as expected from the classical 
micro-turbulent model (Sects. 4 & 5). Thus, we are convinced 
that abundance determination can be possible with the use of the 



" From a detailed measurement of the shifts and shapes of CO lines, 
we showed that the line asymmetries and differential line shifts of CO 
lines are as large as 0.5 km s"' (Paper 111). This results can also be in- 
terpreted as evidence for the presence of the extra molecular layers that 
may have some velocity differences against the photosphere. 

'" We also assume that the MOLsphere is almost transparent to the 
continuum radiation at least in the near infrared and that the infrared 
flux method with the L band flux can still be applied. However, it is 
possible that the MOLsphere will be opaque to the continuum radiation 
at the longer wavelength region as has actually been shown for the radio 
domain of Mira variable stars (Reid & Menten 1997). 
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Fig. 17. Observed abundances of a) carbon (Table 6), b) nitrogen (Aoki & Tsuji 1997; Table 6), c) oxygen (Table 6; If the results 
derived from the H and L band spectra are available, mean values are used), and d) C+N+O, compared with the predictions of 
recent evolutionary models (Claret 2004). The results of FDU and SDU for the case of Z = 0.02 are shown by solid and dash- 
dotted lines, respectively, in each panel. The observed C + N + O abundances are smaller by about 0.16 dex compared with the 
predictions based on Z = 0.02 as can be seen in panel d. This initial metallicity correction of -0.16 dex is applied to the predicted 
C, N, and O abundances in panels a, b, and c, respectively, and the corrected results are shown by dashed lines. In each panel, the 
solar abundances are shown by open square (Anders & Grevesse 1989; Ayres et al. 2006) and by open diamond (Allende Prieto et 
al. 2002). 



weak lines of log Wjv % -4.75. With this caution, our results 
agree well with those by Smith & Lambert (1990), who already 
used only weak lines in their abundance analyses. Then we com- 
pare our results with the predictions of recent evolutionary mod- 
els. 

We plot the resulting C and O abundances in Figs. 17a and 
17c, respectively. Also, we plot N abundances based on our pre- 
vious analysis of the same sample (Aoki & Tsuji 1997), cor- 
rected for the effect of the revised C abundances (14-th column 
of Table 6), in Fig. 17b. For 5 late M giants for which the N abun- 
dances based on the ro-vibrational lines of NH are available, we 
plotted the mean log An values based on the lines of CN and NH. 
Given C, N, and O abundances, we obtain the sum of the C, N, 
and O abundances as plotted in Fig. 17d. The results are com- 
pared with the predictions of the standard evolutionary models 
by Claret (2004): The results of the first dredge-up (FDU) are 
shown by solid lines and those of the second dredge-up (SDU) 
by dash-dotted lines in Figs. 17a-d. 

We first examine the sum of the C, N, and O abundances in 
Fig. 17d. The model prediction simply confirms that the CNO 
processing conserved the sum of the C, N, and O abundances 



corresponding to Z = 0.02, as assumed in the evolutionary mod- 
els. The observed values are generally smaller than the values 
of the models which are log (Ac + An + Aq) ~ 9.16. The mean 
of the observed values in low mass stars (M < 2.5Mq) is log 
(Ac + An + Aq) ~ 9.00 as shown by the dashed line in Fig. 17d, 
and thus the mean metalliciy of the low mass stars we analyzed 
is smaller than the value assumed in the evolutionary models 
by Claret (2004). This metallicity difference should be kept in 
mind in comparing the observed and predicted C, N, and O abun- 
dances. Also, a recent solar value of the C+N+O abundances 
based on the C and O abundances by Alende Prieto et al. (2002) 
is shown by the open diamond and a classical value (e.g. Anders 
Grevesse 1989), recently supported by Ayres et al. (2006), is 
indicated by the open squares in Fig. 17d. The N abundance is 
taken from Anders & Grevesse (1989). 

The observed carbon abundances are generally smaller than 
those of the model predictions, but may be consistent with the 
prediction after the metallicity difference between the evolution- 
ary models and actual stars is corrected as shown by dashed 
line in Fig. 17a. The scatters around the dashed line may be due 
to those of the metallicity in our sample, except for a few ob- 
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Fig. 18. a) The observed '^O/'^O ratios are plotted against stel- 
lar mass by filled circles and by the filled triangles if only lower 
limits can be estimated. Predicted '^O/'^O ratios by FDU and 
SDU are shown by solid and dash-dotted lines, respectively, b) 
The observed '^C/'^C ratios are plotted by filled circles against 
stellar mass. Dashed curve is the predicted '-C/'^C ratios by the 
"Compulsory Deep Mixing (CDM: Eggleton et al. 2008)". 



jects. For example, two stars appear above the predicted FDU 
curve: The object showing the largest deviation is a Cet and the 
other one is A Aqr. It is interesting to notice that these two stars 
also show the largest under abundances of N at M s; 3.6Mq in 
Fig. 17b . 

The nitrogen abundances may roughly be consistent with the 
FDU curve with and/or without metallicity correction in general. 
However, the deviations from the prediction appear to be larger 
in a few objects, and it is as if some extra mixing may be needed 
as in the case of the '"C/'-'C ratios discussed in Sect. 7.4. Most 
lines of both CN and NH used are weak lines that satisfy the 
condition of log W/v ^ -4.75, and N abundances may be free 
from the difficulty encountered in the analysis of the CO first 
overtone bands in Paper I (also Sect. 4.1). 

Finally, the oxygen abundances in Fig. 17c show a similar 
pattern as the C-hN-hO abundances in Fig. 17d, and the scatters 
may represent those of the metallicity in our sample. There is a 
tendency for the oxygen abundances to be larger in higher mass 
objects. The resulting oxygen and carbon abundances correlate 
rather well and the mean Ac/Aq ratio for 19 objects for which C 
and O abundances are determined is 0.37 + 0.10 (std. dev.). 



7.4. Carbon and oxygen isotopic ratios 

The resulting '*0/'^0 ratios (Table7) are plotted against the stel- 
lar mass in Fig. 18a by filled circles or by filled triangles if only 
lower limits could be estimated. For comparison, predicted re- 
sults of FDU and SDU by Claret (2004) are shown by soUd 
and dash-dotted lines, respectively. The rather large variation in 
'^O/'^O of low mass stars is well consistent with the prediction 
of the evolutionary models, confirming the previous analysis by 
Dearborn (1992) using the observed data known at that time. 
In the higher mass stars, some objects have quite low '*0/'^0 
ratios, but the general trend that the '^O/'^O ratios are con- 
fined to lower values can be consistent with the predicted ones. 
Considering the possible uncertainties in the stellar masses, the 
observed and predicted ""O/'^O ratios can be regarded as not 
showing a serious contradiction. 

The observed '^C/'^C ratios (Table 7) are shown by filled 
circles and compared with the predicted ones (Claret 2004) in 
Fig. 18b. The observed '^C/'-'C ratios are now around 10 in con- 
trast to our previous '^C/'^C ratios which were mostly below 
10 (Tsuji 2007), but yet too small compared with '^C/'^'C * 20 
predicted by the standard evolutionary models. This contradic- 
tion has been known for a long time and we simply confirm 
this long-standing puzzle. To resolve this dilemma, extra mix- 
ing of the CN-cycle processed material had to be assumed. For 
example, a deep circulation below the bottom of the convective 
zone, referred to as "cool bottom processing", was proposed ( 
Boothroyd & Sackmann 1999). More recently, it was shown by 
Eggleton et al. (2008) that a deep mixing mechanism should op- 
erate in all the low mass giants due to a molecular weight inver- 
sion resulting from ^He (^He, 2p) "^He reaction in the layers just 
above the H-burning shell. The predicted '^C/'^C ratios by this 
"Compulsory Deep Mixing (CDM)" are shown by the dashed 
line in Fig. 18b. Now, the new prediction appears close to the up- 
per boundary of the observed results, and discrepancies between 
observations and predictions are relaxed if not perfectly. 

7.5. The warm molecular envelope and mass loss from red 
giants 

Given that the presence of the warm molecular envelope or 
MOLsphere is a general phenomenon in normal red giant stars, 
the next problem is to clarify the nature of the MOLsphere in 
detail, and to understand the structure of the entire atmosphere 
extending beyond the photosphere. An initial step towards this 
goal is opened by a detailed MHD simulation of the red giant 
wind by Suzuki (2007), who showed that nearly static regions 
of r « 10^ - 10"* K are formed at several stellar radii above the 
photosphere. It is interesting that the mass loss from red giant 
stars with the wind velocities far smaller than the surface es- 
cape velocity can be understood consistently with his simulation 
in which the stellar winds are effectively accelerated from the 
quasi-static region at several stellar radii above the photosphere. 
The quasi-static region predicted by Suzuki's simulation can 
be identified with the warm molecular envelope that we con- 
firmed in several observational studies discussed in Sect. 7.2; we 
hope that the simulation of Suzuki can be extended to include 
spectral modeling so that more detailed confrontation between 
theory and observations can be achieved. We include in Table 3 
the measured data on the stronger lines together with those on 
the weaker lines used in the abundance analysis, and we hope 
that these data can be interpreted correctly with more sophisti- 
cated models in the future. With such a detailed confrontation 
between theory and observations, we expect to have a consistent 
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solution to the long-standing problem on the origin of mass-loss 
from red giant stars for the first time. 

8. Concluding remarks 

We initiated this study with a hope to clarify the nature of the 
difficulty encountered in our previous analysis of CO (Papers 
I & III ), since we felt that something is still not well under- 
stood in abundance analysis of cool giant stars. At first, we 
had to recognize that the line spectra of cool giant stars are in 
fact highly complicated and may consist of at least three differ- 
ent groups of lines that may have different origins: The weak 
lines of log W/v ^ -4.75, the intermediate-strength lines of 
-4.75 < logW/v ^ -4.4, and the strong lines of log W/v > -4.4. 
Of these three groups, only the weak lines could be understood 
well on the basis of the classical micro-turbulent model of line 
formation. On the other hand, it has been known that the strong 
lines could not be understood as being of the photospheric origin 
alone but should include an extra component originating in the 
outer molecular layers (Paper II). 

The most difficult case was the intermediate-strength lines. 
However, we are finally led to conclude that the intermediate- 
strength lines are essentially the same as the strong lines. This 
means that not only the low excitation strong lines but also 
medium strong lines with LEP as high as 2 eV (Fig. 12) are dis- 
turbed by the contamination originating in the outer molecular 
envelope, which we refer to as MOLsphere for simplicity. This in 
turn means that most of the dominant molecular lines observed 
in the infrared spectra of K - M giant stars should be badly dis- 
turbed by the contamination from the extra component beyond 
the photosphere. At the same time, this finding provides further 
support for the presence of the MOLsphere, on the basis of the 
high resolution spectra of dozens of red giant stars. Thus, the dif- 
ficulty of the abundance analysis we have encountered is nothing 
but a warning that the spectra of cool luminous stars cannot be 
analyzed by the conventional methods using the classical model 
photospheres. Instead, the influence of the extra constituent of 
the atmosphere, which we tentatively referred to as MOLsphere, 
should be quite essential not only in the strong low excitation 
lines but also in many weaker lines that dominate the infrared 
spectra. 

Now we had to know that the largest fault in our previous 
analysis (Paper I) was to have assumed that the critical value 
discriminating the weak and the intermediate-strength lines is 
log W/v X! -4.5 rather than -4.75. For this reason, many lines 
which are disturbed by the contamination from the warm molec- 
ular envelope were included in the abundance analysis of Paper I. 
It should certainly be difficult to attempt abundance analysis on 
such lines including the contamination of the non-photospheric 
origin by the use of a conventional model photosphere, and the 
reason why abundance analysis of cool giant stars often encoun- 
tered difficulties appeared to be quite simple, at least partly. This 
result also implies that it should be difficult in principle to apply 
low resolution infrared spectra to abundance analyses of cool 
luminous stars, and the hybrid nature of the infrared spectra re- 
inforces the reason why high resolution spectroscopy should be 
indispensable in abundance analyses. 

On the other hand, it is fortunate that the weak lines are 
relatively free from disturbance by extra component of non- 
photospheric origin and we are convinced that the abundance 
analysis of cool giant stars can be possible with these lines. With 
a limited number of the weak lines, we could determine the C, 
O, and their isotopic abundances with reasonable accuracy. The 
resulting CNO and their isotopic abundances are roughly consis- 



tent with the predictions of the first and second dredges-up dur- 
ing the red giant phase. We noticed, however, that there are some 
deviations from the evolutionary models especially in '^C, '^C, 
and '"^N abundances in some stars, and it seems that some addi- 
tional episodic events may be needed to explain them. Since the 
CNO isotopes in the universe, especially ^^C, '"^N, and ^^O, are 
largely supplied by red giant stars, further refinements of both 
abundance analyses and evolutionary modelings should be pur- 
sued. 

In conclusion, it turns out that high resolution infrared spec- 
tra include wealth of information not only on photospheric abun- 
dances but also on the structure of the entire atmosphere ex- 
tending beyond photosphere. For this reason, the infrared spectra 
must be of hybrid nature with contributions from multiple com- 
ponents. In the interpretation and analysis of the infrared spectra 
of cool luminous stars, it is essential to keep their hybrid nature 
in mind. 
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Table 2. Observed spectra 



obj. 


Filter 


res.(mK) 


S/N 
60 


Date of obs. 


FTS sp. in ftsdata/ 


alf Tau 


H 


46.19 


Aug. 17, 1976 


kp0222 


alf Tau 


K 


46.19 


76 


Sep. 28, 1976 


kp0226 


alf Tau 


L 


49.90 


159 


Oct. 21, 1977 


kp0628 


del Oph 


H,K 


49.70 


75 


Jun. 24, 1977 


kp0201 


nu Vir 


H 


41.21 


51 


Apr. 08, 1982 


kp0321 


nu Vir 


K 


41.20 


115 


Apr. 09, 1982 


kp0332 


alf Get 


H,K 


46.19 


69 


Sep. 28, 1976 


kp0221 


sig Lib 


H,K 


46.19 


89 


Aug. 18, 1976 


kp0228 


lam Aqr 


H,K 


46.19 


61 


Jul. 06, 1976 


kp0313 


bet Peg 


H,K 


48.76 


103 


May 12, 1976 


kp0104 


tau4 Eri 


H,K 


46.19 


59 


Sep. 29, 1976 


kp0127 


mu Gem 


H 


46.19 


- 


Nov. 07, 1976 


kp0219 


mu Gem 


K 


16.53 


80 


Jan. 27, 1988 


kp0902 


del Vir 


H 


46.19 


54 


Jun. 18, 1976 


kp0212 


del Vir 


H,K 


49.70 


86 


Jun. 24, 1977 


kp0207 


lODra 


H,K 


49.82 


61 


Jun. 28, 1977 


kp0131 


rho Per 


H 


29.70 


44 


Apr. 07, 1982 


kp0317 


rho Per 


K 


22.28 


66 


Apr. 08, 1982 


kp0329 


BS6861 


H,K 


46.19 


42 


Sep. 29, 1976 


kp0114 


del2 Lyr 


H,K 


46.19 


78 


Apr. 11, 1976 


kpOlll 


RRUMi 


H,K 


46.19 


61 


Aug. 18, 1976 


kp0227 


alf Her 


H,K 


49.70 


106 


Jun. 24, 1977 


kp0202 


alf Her 


L 


16.53 


45 


Jul. 09, 1987 


kp0940 


OP Her 


K 


16.53 


142 


Jul. 01, 1988 


kpll24 


XYLyr 


K 


16.53 


86 


Jul. 08, 1987 


kp0930 


RLyr 


H 


29.70 


70 


Apr. 08, 1982 


kp0325 


RLyr 


K 


16.53 


65 


Jul. 07, 1987 


kp0923 


RLyr 


L 


16.53 


107 


Jul. 09, 1987 


kp0941 


RZAri 


H,K 


46.19 


88 


Aug. 17, 1976 


kp0115 


RZAri 


L 


16.53 


65 


Oct. 02, 1987 


kpl028 


30g Her 


H,K 


49.10 


102 


Jun. 30, 1977 


kp0128 


30g Her 


L 


16.53 


- 


Jul. 09, 1987 


kp0939 


SWVir 


H,K 


46.19 


65 


Jun. 18, 1976 


kp0120 


SWVir 


L 


16.53 


88 


Jul. 09, 1987 


kp0936 


RXBoo 


H,K 


48.10 


62 


Jan. 14, 1976 


kp0105 


RXBoo 


L 


16.53 


108 


Jul. 09, 1987 


kp0937 


alfGMa 


H 


29.70 


- 


Apr. 07, 1982 


kp0318 


alfGMa 


K 


22.28 


- 


Apr. 09, 1982 


kp0330 


IRG+ 10216 


L 


16.53 


- 


Oct. 02, 1987 


kpl033 



Table 3. Measured equivalent width data (only a few lines at the beginning of table3 as examples) 



star 


isotope 


v' v" 


rot. tr. 


nu(cm-l) 


loggf 


L.E.R(cm-l) 


logEW/ni 


alf Tau 


C12016 


3 1 


P2 


4199.480 


-6.282 


2154.701 


-4.739 


alf Tau 


C12016 


20 


P9 


4222.953 


-6.128 


172.978 


-4.585 


alf Tau 


C12016 


42 


R68 


4237.805 


-4.241 


12982.145 


-4.721 


alf Tau 


C12016 


3 1 


R8 


4238.289 


-5.600 


2280.402 


-4.611 


alf Tau 


C12016 


42 


R36 


4243.637 


-4.604 


6763.398 


-4.608 


alf Tau 


C12016 


3 1 


R18 


4266.090 


-5.248 


2794.075 


-4.572 


alf Tau 


C12016 


3 1 


R22 


4275.180 


-5.154 


3105.648 


-4.553 


alf Tau 


C12016 


20 


R92 


4292.047 


-4.817 


16004.551 


-5.067 


alf Tau 


C12016 


3 1 


R32 


4292.750 


-4.969 


4148.160 


-4.563 


alf Tau 


C12016 


3 1 


R35 


4296.566 


-4.923 


4533.898 


-4.595 



